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25 BACKGROUND 

The present invention relates to methods for predicting the effects of climate and other 
environmental conditions on the well-being of animals. Specifically, it relates to a method for 
accurately calculating the spatial-temporal effects of a variety of environmental conditions on 
animal individual, population and community dynamics, given the animal's temperature- 
30 dependent behaviors, morphology and physiology, by running integrated microclimate and 
animal models to calculate the discretionary energy and water available to the animal and its 
activity time. 
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Animals are affected by the climate and other environmental conditions in which they 
live. Climate affects animal heat and mass balance and, consequently, body temperature. 
Spatial and temporal variation in the thermal environment provides a template on which 
5 physiological, population and community dynamics are played out. A current challenge in 
biology is understanding how and when those environmental conditions affect developmental 
processes and function during the lives of animals. An ability to calculate those effects would 
enable more accurate prediction and informed decision-making in the management of animal 
environments. Current trends in global climate shifts, land use and changes to other 
10 environmental conditions such as diseases and increases in pesticide use, drive the need for a 
means to accurately predict their effects to prevent, or at least diminish, their negative impacts 
on the well-being of animals. The problem is how to model the spatial -temporal effects of 
climate and other environmental conditions on animals to enable the accurate calculation of 
those effects. 

1 5 Explicit calculations of how climate affects animal heat and mass balance, and the 

consequences for body temperature, were first made in the 1960's (Norris, K.S. 1967. Color 
adaptation in desert reptiles and its thermal relationships. In: Symposium on lizard ecology, 
pp. 162-229. U. of Missouri Press, Columbia, Missouri; Bartlett, P.N. and D.M. Gates, 1967. 
The energy budget of a lizard on a tree trunk. Ecology 48:31 5-322). Those early models 

20 were limited by the lack of models for distributed heat generation internally, distributed 

evaporative water loss internally, and a first principles model of gut and lung function. Also 
missing were a first principles model of porous insulation for fur or feathers, an appendage 
model, and a general microclimate model that could use local macroclimate data to calculate 
the range of local microenvironments above and below ground. 

25 Since then, several developments have occurred to improve our ability to model the 

effects of environmental conditions on animals. These include the estimation of convection 
heat transfer properties from the volume of an animal (Mitchell, J. W. 1976. Heat transfer 
from spheres and other animal forms. Biophys. J. 16:561-569), a countercurrent heat 
exchange model for appendages and modifications based on appendage shapes (Mitchell, J. 

30 W. and G.E. Myers, 1968. An analytical model of the counter-current heat exchange 

phenomena. Biophys. J. 8(8):897-91 1; Wathen, P. et al, 1971. Theoretical and experimental 
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studies of energy exchange from jack rabbit ears and cylindrically shaped appendages. 
Biophys. J. 1 1(12): 1030-1047; Wathen, P. et al, 1974. Heat transfer from animal 
appendage shapes-cylinders, arcs, and cones. Trans, of the ASME. November. 40:536-540), 
the ability to understand the effects of convective heat transport (Kowalski, G J. and J. W. 
5 Mitchell, 1976. Heat transfer from spheres in the naturally turbulent, outdoor environment, 
J. Heat Transfer 98(4):649-653), the development of a general purpose microclimate model 
(Beckman, W. A. et a/. , 1 97 1 . Thermal model for prediction of a desert iguana 's daily and 
seasonal behavior. Trans. ASME Paper No. 71-WA/HT-35:l-7; Porter, W.P. et al, 1973. 
Behavioral implications of mechanistic ecology: Thermal and behavioral modeling of desert 

10 ectotherms and their microenvironment. Oecologia 13: 1-54; Mitchell, J.W. et al (eds) 1975, 
Microclimatic modeling of the desert, pp. 275-286. Scripta Book Co., Washington, D.C.), the 
ability to calculate the percent of thermally available habitat (Grant, B.W. and W.P. Porter 
1992. Modeling global macroclimatic constraints on ectotherm energy budgets. Am. Zool. 
32:154-178), the development of porous media models (Kowalski, G. J. 1978. An analytical 

1 5 and experimental investigation of the heat loss through animal fur. In: Department of 

Mechanical Engineering. University of Wisconsin, Madison, Wisconsin.), the extension of 
the models to radial instead of Cartesian coordinates and the implementation of first principles 
fluid mechanics in the porous media (Stewart, W.E. etal, 1993. Prediction of forced 
ventilation in animal fur under ideal pressure distribution. Functional Ecology 7:487-492; 

20 Budaraju, S. et al , 1994. Prediction of forced ventilation in animal fur from a measured 

pressure distribution. Proc. Roy. Soc. London B 256:41-36; Budaraju, S. et al, 1997. Mixed 
convective heat and moisture transfer from a horizontal furry cylinder in transverse flow. Int. 
J. Heat and Mass Transfer 40:2273-2281), and a test of an ecototherm and microclimate 
model to estimate a species' survivorship, growth and reproduction at a continental scale 

25 (Adolph, S.C. and W.P. Porter, 1993. Temperature, activity, and lizard life histories. Am. 
Nat. 142:273-295; Adolph, S.C. and W.P. Porter, 1996. Growth, seasonality, and lizard life 
histories: Age and size at maturity. Oikos 77:267-278). However, none of these attempts 
have accomplished the goal of creating a fully integrated set of models that incorporate all of 
the conditions needed to accurately predict how animals (both ecototherm and endotherm) 

30 will react to changes in their environment. 
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SUMMARY 

The present invention is directed to overcoming these limitations by providing a 
system and method to enable users to accurately calculate the effects of climate and other 
environmental conditions, and the spatial and temporal variation thereof, on the temperature 
dependent behavior and physiology of animals, and the subsequent effects of those same 
conditions on population dynamics and community structure. 

Accordingly a methodology for accurately predicting the spatial-temporal effects of 
environmental conditions on animals is presented. Several objects and advantages of the 
present invention are: 

a) to provide improved and fully integrated microclimate and animal models which may be 
used in coordination to accurately predict the effects of climate and other environmental 
conditions (including artificial environments) on animals, requiring relatively few, 
easily measured data to estimate food and water requirements, potential for activity 
time, growth, and reproduction for a wide variety of habitats; 

b) to provide a methodology that can generate the information necessary to identify 
potential reserves and/or transplantation sites, and to modify and/or manage existing 
habitats to optimize the growth, reproduction and survivorship of threatened or 
endangered species; 

c) to provide a method that enables the incorporation of the spatial-temporal variation in 
vegetation and topography in the analysis of the effects of environmental conditions, 
such as changes in global climate, land use, diseases and toxicants, on animals; 

d) to provide a means to accurately calculate discretionary energy and water from first 
principles, thus enabling the calculation of growth and reproduction under any kind of 
environmental circumstance, and under any kinds of animal characteristics; 

e) to provide a microclimate model with a solar radiation model incorporated therein, the 
solar radiation model enabling the input of solar radiation information from the 
geographic coordinates of the locality, thus eliminating the need for direct 
measurements at the site, the solar radiation model further allowing for the calculation 
of the effects of slope on the amount of solar radiation incident on the surface and for 
the calculation of solar radiation when the sun is positioned below the local horizon; 
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f) to provide a microclimate model able to calculate the effects of vegetation on the 
temperature and wind profiles, topography on air temperature profiles over elevation, 
and the effects of slope on the temperature; 

g) to provide an interactive combination of heat and mass transfer engineering and specific 
aspects of an animal's morphology, physiology and temperature-dependent behaviors, 
to calculate an animal's preferred activity time (i.e., the time that minimizes size- 
specific heat/water loss); 

h) to provide an animal model that incorporates a porous media model for insulation 
improved by the addition of a solar radiation component enabling the user to calculate 
the absorbed solar radiation when an animal is in the sunlight; 

i) to provide a porous media model further improved by the addition of an appendage 
model for birds and other animals to account for energy fluxes at their appendages; 

j) to provide an animal model further improved by the incorporation therein of molar 
balance models for the lungs and gut to allow for the coupling of food ingested to 
respiratory gas exchanges and evaporative water loss from the respiratory system; 

k) to provide a means of estimating the food intake of animals in natural circumstances 
using a gut model based on first principles rather than regression analyses; 

1) to provide an improved means of calculating heat loss from an endotherm by 
incorporating the effects of solar radiation and relative humidity as factors in the 
calculation; 

m) to provide an animal model further improved by the incorporation therein of a 

thermoregulatory model that yields metabolic calculations as a function of temperature 

and incorporates temperature regulation responses by the animal to maintain a core-skin 

temperature gradient above a user-specified minimum value; 
n) to provide a means of implementing the models such that a user may run the models 

with access to large databases of environmental and animal data on a variety of 

platforms, including a personal computer; and, 
o) to provide a means by which a user is able to specify the conditions set used by the 

model(s) and/or the results calculated by the model(s), in order to customize the model 

to the user's particular needs. 
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Further objects and advantages of the present invention will become apparent from a 
consideration of the drawings and ensuing description. 

BRIEF DESCRIPTION OF DRAWINGS 

5 In the accompanying drawings: 

FIG. 1, is a flow diagram of the general method for the present invention that also 
represents a logic flow that can be implemented in discrete circuits. 

FIG. 2, is a summary of the general method for the animal model. 
FIG. 3a, is a diagram of a cross-section of an animal with porous insulation and heat 
1 0 fluxes including uniform heat generation per unit volume and uniform heat dissipation by 
respiration per unit volume. 

FIG. 3b, is an illustration of the coupled heat and mass transfer equations and how the 
models may be used to calculate the energetic consequences of climate-disease-toxicant 
interactions and to assess energetic constraints on community structure. 
1 5 FIG. 4, is a representation of the molar balance models of respiratory and digestive 

systems coupled to each other and to metabolic heat demand to maintain core temperature. 
FIGS 5a to 5d, present a key to the variables used in the drawings and algorithms. 
FIG. 6, is a graph showing results generated by the present invention, that show the 
variation in metabolic rate as a function of mass and temperature for the Orange-bellied 
20 Parrot. 

FIG. 7, is a plot of the spatial variation in soil degree days in relation to the range of 
the Chuckwalla lizard, generated by the present invention to discover the effects of 
temperature on the lizard's distribution. 

FIG. 8, is a flow diagram depicting the general system of the present invention as 
25 implemented in a client-server environment. 

FIG. 9, is a spatial plot summarizing the discretionary energy of a 60 kg Mountain 
Lion for the month of July, as a function of slope and aspect of the topography of a portion of 
its range. 

FIG. 10, is a depiction of the basic components of the system of the present invention. 

30 
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OVERVIEW 

The present invention is a method and system for calculating the spatial-temporal 
effects of a variety of environmental conditions (including macroclimate, topography, 
vegetation and the like) on animal individual, population and community dynamics, given the 
5 animal's temperature-dependent behaviors, morphology and physiology. The method for the 
present invention comprises using two integrated models, a microclimate and animal model, 
to translate the environmental and animal conditions into their effects on an animal's 
energetics and behavior. The system of the present invention sets forth a way in which the 
methodology may be employed by users via a wide range of computerized platforms. 

10 Referring to FIG. 1, an overview of the method for the present invention is presented. 

The two integrated models of the present invention include a microclimate model 300 and an 
animal model 500. The microclimate model 300 translates the environmental conditions into 
the microclimatic conditions above and below ground that may be experienced by an animal 
at a specified location. The animal model 500 calculates the effects of the microclimatic 

15 conditions on an animal's metabolic and other variables. Those results are then available for 
further calculations aimed at assessing the ultimate effects of the environmental conditions on 
an animal's physiological performance, geographic distribution and community dynamics. 

The microclimate model 300 of the present invention improves upon a basic model 
previously disclosed in Porter et al 1973 (ref. above) and by McCullough and Porter (1971. 

20 Computing clear day solar radiation spectra for the terrestrial ecological environment. 
Ecology 52(6): 1008-1015). The improvements to the microclimate model 300 include the 
incorporation therein of a solar radiation model. The solar radiation model enables the input 
of solar radiation information from the geographic coordinates of the animal's locality, 
eliminating the need for direct measurements at the site. The solar radiation model also 

25 allows for the calculation of solar radiation when the sun is positioned below the horizon in 
addition to when the sun is above. The microclimate model 300 is further improved in being 
able to calculate the effects of vegetation on the temperature and wind profiles, and the effects 
of slope on the temperature, as well as enabling a user to select the heights and depths above 
and below ground at which the microclimates are characterized. 

30 The animal model 500 of the present invention builds upon the fact that in biological 

systems, mass and energy are interconnected. For example, the amount of mass consumed by 
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an animal is related to the amount of chemical energy that can be liberated by metabolic 
processes in the form of heat. Earlier attempts to develop coupled equations of heat and mass 
balance to quantitatively measure the metabolic variables involved in each were reported in 
two earlier works (Porter, W.P. and P. A. McClure, 1984. Climate effects on growth and 
5 reproduction potential in Sigmodon hispidus and Peromyscus maniculatus. In: 4 Winter 

Ecology of Small Mammals.' J. Merritt (ed), at pp. 173-181; W. P. Porter, 1989. New animal 
models and experiments for calculating growth potential at different elevations. 
Physiological Zoology 62: 286-313). These earlier developments illustrate that relatively few 
variables can govern processes at multiple levels of biological organization. For example, 

10 accurate solutions for the discretionary energy and mass of an animal may be used to calculate 
the life history variables of growth and reproduction, and solutions for the total annual activity 
time of an animal may be used to calculate survivorship. The fact that accurate calculations 
of these few metabolic variables can predict effects on animals, underpins the robustness and 
broad applicability of the animal model 500 of the present invention. 

1 5 An earlier attempt to couple an animal model to a microclimate model in order to 

predict the effect of climate and other environmental conditions on animals, was reported in 
W.P. Porter et al (1994. Endotherm energetics: from a scalable individual-based model to 
ecological applications. Aust. J. Zool., 42:125-162 ). The animal model 500 of the present 
invention builds upon this earlier attempt and incorporates several improvements that enable 

20 the model to more accurately predict the effects of environmental conditions on animals when 
integrated with the microclimate model 300, and that enable its application to both 
endothermic and ecotothermic animals. 

These improvements of the animal model 500 of the present invention, are disclosed 
herein and include several additions or modifications to some of the components of the 

25 previous version, including: (1) the incorporation of a solar radiation component into the 

porous media insulation model to enable the calculation of absorbed solar radiation when the 
animal is exposed to light; (2) the inclusion of an appendage model to improve the accuracy 
of the heat balance analysis by inclusion of the heat exchanges at the appendages of large 
birds and mammals; (3) the inclusion of a novel thermoregulatory model that yields metabolic 

30 calculations as a function of temperature and incorporates temperature regulation responses 
(behavioral and physiological) by the animal to maintain a minimum core-skin temperature 
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gradient; (4) the ability to track the time an animal is active/feeding to calculate a total annual 
activity time; and, (5) the inclusion of molar balance models for the lungs and gut to allow the 
coupling of food ingested to respiratory gas exchanges and evaporative water loss from the 
respiratory system. 

5 The animal model 500 calculates heat (i.e., energy) and water mass balances for 

required animal expenditures that are contingent on their selected or maintained core 
temperature and, from that, further calculates the required wet food intake that is constrained 
to be no greater than the body mass per day. Activity time calculations depend on light or no 
light availability and local microclimates that may force heat or water stress and reduce or 
10 limit activity. Together, the subroutines of the model 500 comprise an integrated basic 
principles model for understanding animal energetics. When coupled with Geographical 
Information Systems (GIS) based information on climate, topography and vegetation, the 
integrated microclimate 300 and animal 500 models can be used to generate spatially explicit 
^ interpretations of animal energetics and behavior. 

C 1 5 The method for implementing the present invention provides a means by which a user 

rQ may run the models of the present invention with access to large environmental conditions 
\T and animal databases on a wide variety of computerized platforms, including a personal 
P computer. 

□ While the methods disclosed herein may be described and shown with reference to 

^20 particular steps performed in a particular order, it will be understood that these steps may be 

fy 

y combined, sub-divided, or re-ordered to form an equivalent method without departing from 
r** the teachings of the present invention. Accordingly, unless specifically indicated herein, the 
order and grouping of the steps is not a limitation of the present invention. 

Referring now specifically to the figures, in which identical or similar parts and/or 
25 steps are designated by the same reference numerals throughout, detailed descriptions of the 
method and system of the present invention are given. It should be understood that the 
following detailed descriptions relate to the best presently known embodiments of the 
invention. However the descriptions are not intended to limit the invention to those 
applications. Rather, the present invention is intended to apply to other types of animal 
30 individual, population, and community research, and any other application to which the 

technology may be beneficially applied. For example, the invention can be applied to other 
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conservation applications and to the prediction of global and other climate changes on animal 
distribution, growth and survivorship, in order to develop intervention procedures to ensure 
their well being. The present invention can, therefore, assume numerous other embodiments, 
as will become apparent to those skilled in the art, without departing from the appended 
5 claims. 

DETAILED DESCRIPTION - METHOD OF CALCULATING EFFECTS 

The models of the present invention comprise several algorithms with terms that are 
defined in the following detailed description. Please refer to FIGS. 5a to 5d for a list of terms 
1 0 and their definitions. 

a. Data Inputs 

Referring to FIG. 1 (a flow diagram of the method for the present invention), two sets 
of data are distinguished as input to the microclimate 300 and animal 500 models of the 
present invention. One is a set of environmental conditions data 120 and the other is a set of 
15 animal characteristics data 170. 

i. Set of environmental conditions data 120 
The environmental conditions data set 120 is a collection of data from an 
environmental conditions database 100 and from a user-defined set of variables and 
conditions 150. 

20 The environmental conditions database 100 is comprised of data relating to climate, 

topography, soil type and vegetation per locality and time. These data may be available from 
public sources (such as the National Oceanic and Aeronautics Administration, NO A A) and/or 
other sources [e.g. soil properties and vegetation type data - may be generated by using 
programs designed to do so, from data collected at local biogeography labs, or from other 

25 sources (e.g. Ramankutty and Foley. 1999. Estimating historical changes in global land 

cover: Croplands from 1700 to 1992. Global Biogeochemical Cycles 13:997-1027) ; certain 
climate information - from individual climate stations and/or from global climate data sets 
(such as that of Leemans and Cramer. 1991 . The IIASA database for mean monthly values of 
temperature, precipitation and cloudiness on a global terrestrial grid. Laxenburg, Austria, 

30 International Institute for Applied Systems Analysis); and, spatial locations data - from the 
United States Geological Survey (USGS)]. Climate data may include temperature, rainfall, 
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wind velocity and direction, and solar and infrared (IR) radiation. Vegetation data may 
include type and abundance of plant cover. Topographic data may include the slope and 
aspect of the terrain and soil type. Other types of data may also be input to the models if 
desired by the user and as they become available. 
5 The user-defined variables and conditions 150, input to the environmental conditions 

data set 120, may include a set of localities and definition of time intervals over which the 
models are to be run, the types of vegetation or other food type(s) that comprise an animal 
species' diet, and the preferred temperature range and average height above ground of the 
animal 

10 Localities may be defined by latitude and longitude (which, in 0.5 degree increments, 

define approximately 30 x 30 mi areas), by use of the universal transverse Mercatur (UTM) 
grid system to define smaller areas in the range of 30 x 30 m , or by other means. Other 
locality size increments are also possible. By defining the localities set, a user may determine 
O the geographic scale of the microclimate conditions data 400 and of the analysis generally. 
JU5 For example, the analysis of the present invention may span an animal species' home range at 
m the individual level, the distribution of individuals at a particular population's level, or it may 

span the entire species' distribution worldwide. 
15 Time intervals may likewise be defined so that the models (both microclimate 300 and 

s _ animal 500) compute outcomes over customizable time intervals. For example, the time 
[t]20 intervals may be specified to calculate hourly, daily, monthly, yearly or other time periods 
j ^ according to the user's desired analysis. In one embodiment, the time intervals will be set to 
O analyze conditions during an average day of each month, over each month in a year to analyze 
animal conditions over an entire year. 

ii. Set of animal characteristics data 170 
25 The animal characteristics data set 170 is a collection of data from an animal 

characteristics database 190 and from a user-defined set of variables and conditions 150. 

The animal characteristics database 190, is a compilation of data over many species of 
animals. The animal characteristics database 190 includes data on a species' morphology, 
body allometry (i.e., dimensions, weight, geometry and the like), physiology (including core 
30 temperature) and other information. Morphological data include those relating to spatial 

reflectance and include mean length of fibers (i.e. fur or feathers), fiber density as a function 
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of depth, fiber diameter, depth of insulation, and solar reflectivity and transmissivity of fibers. 

In addition to the data available from the animal characteristics database 190, the user 
may define a set of variables and conditions 150 for the particular species of concern that also 
make up the animal characteristics data set 170. User defined variables and conditions 150 
5 may include the average height of the particular animal species specified, the minimum core- 
skin temperature gradient, the animal's null condition [i.e. when it is active (e.g. whether it is 
diurnal, nocturnal, crepuscular), and/or where it is active (e.g. whether it is arboreal or 
fossorial)], the animal's food type and quality (i.e. the proportions of carbohydrate:protein: 
lipids; percent dry matter), percent digestive efficiency, percent water in feces, percent urea in 

10 urine, and a multiplier from 1 to 7 above basal metabolic energy requirements to specify food 
intake to meet current needs, among others. The user may define whether and what pollutants 
(such as pesticides) are present in the localities, and whether and what other species (e.g., 
competitors, predators, prey) are present. The user, by inputting particular types and values of 
the data, may generate model outputs customized to their particular inquiries. 

1 5 b. Microclimate model 300 

The microclimate model 300 of the present invention functions to translate the set of 
environmental conditions data 120 (including large-scale macroclimatic information) into a 
description of the microclimates experienced by the animal The model 300 is based on a 
basic model (hereinafter referred to as "basic model") previously disclosed as described 

20 below, but further incorporates several improvements thereto. 

The improvements to the microclimate model 300 include the incorporation of a solar 
radiation model. The solar radiation model enables the input of solar radiation information 
from the geographic coordinates of a locality, eliminating the need for direct measurements at 
the site. The solar radiation model also allows for the calculation of solar radiation when the 

25 sun is positioned below the horizon. The microclimate model 300 is further improved in 

being able to calculate the effects of vegetation on the temperature and wind profiles, and the 
effects of slope on the temperature, as well as enabling a user to select the heights and depths 
above and below ground at which the microclimates are characterized. 

The basic model calculates certain key climate variables from basic principles. The 

30 basic model is a one-dimensional finite difference model that uses data from the set of 
environmental conditions data 120 (including air temperature and wind speeds at certain 
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heights, average annual temperature to specify deep soil temperature, soil properties, humidity 
at a known height, measurements of solar radiation, and geographic coordinates) to 
simultaneously solve the heat and mass balance equations for the ground surface and below 
(solving the equation numerically with an Adams Predictor-Corrector algorithm), and thereby 
5 the above- and below-ground microclimates, over specified time intervals (e.g., hourly 

simulations for the average day for each month of the year for any locality). The basic model 
also computes wind speed and temperature profiles from the ground surface (for flat surfaces 
only) to a two meter reference height (i.e. the height where meteorological data are typically 
measured). Clear sky solar radiation (when the sun is above the horizon) is measured and 
1 0 these values input to the basic model. 

Certain of the basic model calculations were disclosed in prior publications. Air 
temperature and wind speed profile, soil temperature, and soil surface and below boundary 
conditions calculations, are set forth in Porter et al [1973. Behavioral implications of 
O mechanistic ecology: thermal and behavioral modeling of desert ectotherms and their 
/Si 15 microenvironment. Oecologia 13:1-54. Wind velocity calculation set forth at equation (2) on 
^ p.4; air temperature calculation set forth at equation (3) on p. 5; soil temperature calculation 
y* set forth at equation (4) on p. 5; soil surface and below boundary conditions are set forth at 

equations (5) to (1 1) on p. 6], the disclosure of which is incorporated herein by reference. 
- Available solar radiation calculations when the sun is positioned above the horizon are 

gfiO described in McCullough and Porter [1 97 1 . Computing clear day solar radiation spectra for 
v% the terrestrial ecological environment. Ecology 52(6):1008-1015. Solar radiation is 
O calculated according to the flow chart (Fig. 2) on p. 1013, which refers to the methodology 
and the various equations used in the calculations employed], the disclosure of which is 
incorporated herein by reference. 
25 The basic model uses these data, together with information particular to the animal 

species selected (i.e. average height above ground and its preferred temperature range; both of 
which are specified in the set of animal characteristics data 190), to calculate the climate 
boundary conditions for an average individual of the species (including solar radiation, IR 
fluxes from the sky and ground, air temperature, wind speed, and relative humidity of the air 
30 passing over the animal) and the percentage of habitat thermally available to it. The 

percentage of thermally available habitat is the percentage that allows the animal to stay 
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within its preferred temperature range, thereby avoiding significant thermoregulatory heat 
stress costs. 

As mentioned above, the microclimate model 300 of the present invention provides an 
improved version of the basic model including: (1) the incorporation of a solar radiation 

5 model; (2) the ability to calculate the effects of vegetation and slope on temperature and wind 
profiles; (3) the ability to calculate the effects of topography on air temperature profiles over 
elevation; and, (4) the ability for a user to specify the heights and depths above and below 
ground at which the microclimates are characterized. 

L The solar radiation model. 

10 The microclimate model 300 improves upon the basic model in that it incorporates 

therein a solar radiation model to calculate the amount of solar radiation incident on the 
ground. The solar radiation model of the present invention calculates the amount of solar 
radiation absorbed per unit area of ground, duration of solar radiation on the ground, and 
spectral quality, if needed. The basic model, by contrast, takes input on solar radiation 

1 5 measured and/or calculated elsewhere. 

Data used by the solar radiation model are provided by a database of solar radiation 
defined by spectral bandwidth energies available above the earth's atmosphere. Current Julian 
date is used to compute earth-sun distance to adjust solar energy per unit area. Solar radiation 
at the earth's surface is calculated by the extinction properties of the atmosphere, depth of the 

20 atmosphere and geographic coordinates of a locality at a given time. This eliminates the need 
for direct measurements at a site. 

The solar radiation model is also improved in that it accounts for the effects of slope 
on solar radiation reaching the ground, whereas the former solar radiation calculations 
assumed flat ground surfaces only. The effects of slope on the amount of solar radiation 

25 reaching a surface may be significant. Accounting for those effects increases the model's 
accuracy in solar radiation calculations, and thereby also calculations based upon them, such 
as ground and air temperatures. 

A further improvement of the solar radiation model is that it is now able to calculate 
the amount of solar radiation when the sun is below the horizon, enabling more accurate 

30 calculations for the amount of solar radiation present when crepuscular animals are active at 
dusk and dawn. 
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The solar radiation model and its various improvements are incorporated into the 
microclimate model 300 as follows. 

Solar radiation dataset. Solar radiation data used by the solar radiation model 
component of the microclimate model 300, are calculated using a database of solar radiation 
spectral bandwidth energies available above the earth's atmosphere (F.S. Johnson (private 
communication to inventor, in Ch.l, p. 2 of Robinson, N. 1966. Solar Radiation. Elsevier Publ. 
Co., New York. 347p.). Current Julian date is used to compute earth-sun distance to adjust 
solar energy per unit area on top of the atmosphere. Solar radiation at the earth's surface is 
calculated from the extinction properties of the atmosphere {see, McCullough and Porter, 
1971, ref. above), depth of the atmosphere and geographic coordinates of a locality and time. 
Measurements of solar radiation taken at the actual location are no longer required, although 
they may be used for verification purposes. The database of spectral bandwidth dependent 
atmospheric extinction coefficients and solar radiation spectral bandwidth energies is 
contained in a data file, Solar.dat (not shown). 

Effect of slope on solar radiation. The solar radiation model further improves the 
basic model by incorporating the effects of slope on the amount of solar radiation reaching the 
ground. The effects of slope are calculated by equations based on work by Hosmer, G.L. 
(1910. Practical Astronomy. Wiley & Sons. New York, or Hosmer, G.L. and J. M. Robbins. 
1948. Practical Astronomy. Wiley & Sons. 355 pp.). The summary equation is presented in 
W.D.Sellers (1965. Physical Climatology, University of Chicago Press, 272 pp. Equation 
No. 3.15 on p. 35), the disclosure of which is incorporated herein by reference). 

Below horizon solar radiation. The solar radiation model enables the calculation of 
the amount of solar radiation when the sun is below the horizon, as at dusk and dawn. Prior to 
this the basic model assumed that the sun was full sun and located above the horizon. Values 
for below horizon solar radiation enable more accurate measures for the amount of solar 
radiation in the microclimatic environments of crepuscular animals. 

The solar radiation model of the present invention calculates the amount of solar 
radiation when the sun is below the horizon by first computing the amount of skylight {see, 
G.V. Rosenberg, 1966. Twilight. Plenum Press, pp. 18-19). First, the total lumens are 
calculated from a regression of data from "Twilight" (which is Elog = 41.34615384 - 
0.423076923 * Zenith angle). The sky luminosity is equal to 10**Elog. Lumens/m2 are then 
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converted to W/m2 (see, Documenta Geigy Scientific Tables. 1966. 6th edition. K. Diem, 
ed., p. 239), by multiplying by the mechanical/electrical equivalent of light (i.e., 1.46*10 " 3 
kW per lumen). 

ii. Effects of vegetation on microclimate. 
Vegetation-climate interactions can impact directly and indirectly animal food and 
water requirements, growth and reproduction potential, activity patterns, population 
dynamics, and structure of food webs. For example, loss of vegetation in cold climates can be 
particularly challenging for large animals that depend on intact evergreen vegetation for 
radiant protection during the cold winter nights that they must endure. The effects of 
vegetation on an animal's microclimate are calculated by the microclimate model 300 of the 
present invention. 

The microclimate model 300 calculates the effects of vegetation on the temperature 
and wind profiles from the ground level to a 2 m reference height, including at the animal's 
average height above ground. 

Under dense vegetation, the overhead sky radiant temperature (Tsky) is assumed to be 
equal to the local 2 m air temperature (Tair,2m). The leaves or needles of trees and/or bushes 
may provide radiant insulation because they generally have high heat transfer coefficients, 
which causes them to be convectively warmed at night or cooled during the daytime. 
Multiple layers of vegetation above 2 m may provide further radiant insulation. 

Modification of wind velocity profiles by vegetation may be measured empirically. In 
the absence of such information, the "bushel basket" formula can be used to estimate velocity 
profiles between objects of known dimensions and distances (C.R. Stearns. 1967. 
Micrometeorological studies on the coastal desert of Peru. Ph.D. thesis, University of 
Wisconsin, Madison, WI). 

Whether temperatures between the 2 meter height (Tair,2m) and ground surface 
(Tgrd) vary, depends upon the amount of sunlight striking the ground. One extreme is 
complete shade, which results in no temperature differences between ground and two meters. 
The other extreme is no shade, which results in a typical logarithmic air temperature profile. 
Changes in the percent of ground shade alter ground surface temperature, which in turn, alters 
the temperature profile between the surface and 2m shade temperature. Though vegetation is 
assumed to alter wind velocity profiles, temperature profiles are currently assumed to be 
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unaffected by vegetation, except for the degree of ground shade specified by the user. 

iii. Effects of slope on microclimate. 

The amount of solar radiation reaching the ground surface depends on the slope of the 
ground. Hence an improvement over the prior basic model, is the incorporation of slope into 
5 the present microclimate model 300 to improve the accuracy of solar radiation calculations. 
This is done by equations based on work by G.L. Hosmer (the summary equation is No. 3.15 
on p. 35 of Physical Climatology by W. D. Sellers. 1965. University of Chicago press. 272 
pp., and is incorporated herein by reference). 

iv. Effects of topography on microclimate. 

1 0 Topography can affect air temperature profiles over elevation. The present 

microclimate model 300 is able to calculate those effects by using adiabatic principles to 
adjust temperatures based on elevation differences from temperatures at a known height. 
There are various corrections used from 9.9 degrees C. per kilometer in non-mountainous 
regions, to 5.5 degrees C. per kilometer in mountainous regions. The equation using meters 
15 and degrees C. is: 

current temperature = reference temperature at reference elevation + 0.0055 * 
(reference elevation - current elevation). 

v. User-specifiable heights and depths. 

The microclimate model 300 has also been improved to allow a user to specify the 
20 heights and depths above and below ground at which the microclimates are characterized. 
The user is thus able to customize the microclimate outputs to increase the accuracy for a 
particular species, to enable a broader array of scenarios to be posed in order to better manage 
an animal population, or to achieve any number of other goals. From a knowledge of the 
animal morphology, the user specifies the animal's average height above ground. The animal 
25 height is then used by the program to compute the wind speed and temperature for profiles 
that the user has specified, by selection of a reference height for the null calculations for an 
animal at each locality. 

vi. Microclimate model conditions output 400 

The microclimate model conditions output 400 includes the solutions to the heat and 
30 mass balance equations for the ground surface and below, the climate boundary conditions 
experienced by an individual of the animal species selected, and the percentage of habitat 
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thermally available to it. 

Many of the climate boundary conditions calculated by the microclimate model 300 
are illustrated in FIG. 3a, a diagram of a cross-section of an animal with porous insulation 
and heat fluxes including uniform heat generation and uniform heat dissipation by respiration. 
5 They include the clear sky solar radiation (Qsolar), IR fluxes from the sky (QIR,sky) and 
ground (QIR,grd), air temperature at the average height of the animal (Tair), ground 
temperature (Tgrd), and the wind speed (V) and relative humidity (RH) of air passing over the 
animal at its average height above ground. 

The percentage of habitat thermally available to the animal is also calculated based on 
10 the climate boundary conditions and the specified preferred temperature range of the animal. 

c. Animal model 500 

The animal model 500 of the present invention uses the microclimate conditions 
output 400 and data included in the set of animal characteristics 170, to calculate an animal 
conditions output 600 (see, W.P. Porter et al 2000. Calculating climate effects on birds and 
y315 mammals: impacts on biodiversity, conservation, population parameters, and global 
g community structure, Amer. Zool. 40:597-630). In general, the output 600 includes the values 
for several energetic variables including total discretionary energy and water, the temperature- 
f?| dependent activity of the animal and the total annual activity time of the animal. Using the 
^ values for these energetic variables, spatial analyses and further calculations 800 may be 
yy 20 performed to analyze spatial variation and to solve for a variety of population and community 
y level variables such as survivorship, growth and reproduction. 

fT In overview, the animal model 500 calculates the aforementioned energetic variables 

by performing coupled analyses of heat and mass balance and by calculating the total annual 
activity time of the animal. Different subroutines of the animal model 500 are employed to 

25 perform various aspects of the heat and mass balance analyses and to track activity time of the 
animal The subroutines include a porous media model for insulation, an optional appendage 
model, coupled models of the gut and lungs, and a thermoregulatory model. 

Many of the subroutines are run in parallel and/or iteratively to solve for the various 
values that drive the final solutions. Together they enable the animal model 500 to accurately 

30 calculate the amounts of discretionary energy (Qdiscr) and water (Mw,discr) from first 
principles (i.e., rather than regression analyses), thus enabling the accurate calculation of 
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growth and reproduction under any kind of environmental circumstance, and under any kind 
of animal property. 

Referring to FIG. 2, the animal model 500 portion of FIG, 1 is expanded into a flow 
chart summary of its basic processes. The following description explains how the various 
5 sub-models and other calculations performed by the animal model 500, contribute to each of 
the steps indicated in FIG. 2. 

The heat and mass balance analyses performed by the animal model 500 represent the 
coupled mechanisms of heat and mass transfer. Referring to FIG 3b, coupled heat and mass 
transfer equations are represented that illustrate how they are used to quantitatively evaluate 
10 energetic consequences of environmental interactions (e.g., with climate, disease, toxicants 
and the like), and to assess energetic constraints on community structure. 

The equations shown in FIG 3b share common terms in the heat and mass transfer 
equations for metabolism and evaporation. The equations emphasize that the discretionary 
^ mass (Mdiscr) available for growth (Mgr), reproduction (Mrep), and fat storage (Ms), depends 
yQ 15 upon (1) environmental conditions that affect heat fluxes in the diagonal heat balance 
^ equation, (2) mass ingested and absorbed which determines the amount of dry matter and 
V water ingested and absorbed, and (3) metabolic rates and evaporative water losses. Metabolic 
fll rates and evaporative water losses are affected by an animal's physiology (e.g.., body 
L. temperature maintained) and behavior (e.g., microclimates selected, activity levels, locomotor 
1320 costs to obtain food with distribution that changes, appetite levels and the like). 
Lyj Heat and mass transfer may occur through the flesh, fat and porous insulation layer 

J=f (i. e., fur or feather) of animals. Referring to FIG. 3a, a diagram of a cross-section of an 
animal with a porous insulation layer and heat fluxes including uniform heat generation 
(Qgen) and uniform heat dissipation by respiration (Qresp;) is depicted [W.P. Porter et al, 
25 2000. (ref. above), Fig. 2 on p. 602]. The cross-section could represent a torso whose 

geometry may be approximated by a cylinder, sphere, or ellipsoid, or even a cross section 
through an appendage, if the heat loss by respiration is removed (see appendage model 
below). There may or may not be porous insulation beyond the skin. 

FIG. 3a depicts the information needed for the heat and mass balance analyses 
30 performed by the animal model 500. Data needed are provided in the set of animal 

characteristics data 170 [including, mean length of the fibers (hair or hair-like elements in 
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feathers); fiber density as a function of depth, fiber diameter, the depth of the insulation, and, 
if the animal is diurnal and exposed to sunlight, solar reflectivity and transmissivity of the 
fibers]. Also needed are data from the microclimate conditions 400 that specify the climate 
boundary conditions for an individual [solar radiation (Qsolar), IR fluxes from the sky 
5 (QIR,sky) and ground (QIR,grd), air temperature at the animal's average height (Tair), ground 
or other surface temperature (Tgrd), and the wind speed (V) and relative humidity (RH) of air 
passing over the animal at its average height above ground]. The climate boundary conditions 
are calculated by the microclimate model 300 based on the animal's average height above 
ground and the microclimate calculations for environmental conditions above ground (see 
10 section on microclimate model 300 above). 

i. Performing the heat balance analysis 501 
Performing the heat balance analysis 501 comprises running a porous insulation model 
to simultaneously calculate the amount of heat entering the animal from external sources 
u (Qin), and the amount of heat leaving the animal to the environment (Qout), contingent on the 
yyl 5 current guess for the metabolic heat generation needed to maintain the animal's core 
^ temperature (Qgen). Simultaneous calculations are required because in porous media, diffuse 
H radiation fluxes and conduction interact with each other to affect the insulation temperature 
p profile from skin to insulation-air interface. If the core-skin temperature gradient is too small, 
L. then a thermoregulatory model 537 is invoked until a stable solution is found that will allow a 
CQ20 core-skin temperature gradient of at least 0.5 degrees C to exist. 

lj 1. Calculating the amount of solar radiation absorbed by the animal 

3 (Qsolar,a) 510. 

Computing the amount of solar energy absorbed (Qsolar,a), is accomplished as 
follows. 

25 Direct beam solar radiation (Qsolar,dir): 

Qsolar,dir = ABSAN * ASILN * (LOO - PCTDIF) * Qnorm * (100 - ashade) / 
100 
where, 

ABSAN is animal absorptivity, 
30 ASILN is silhouette area normal to the sun, 

PCTDIF is percent diffuse solar radiation, 

(Express Mail label # EK796645358US, Filing Date 23 AUG 2001)2Q 



Qnorm is total solar radiation incident on a surface normal to the sun's 

direct beam, and 
ashade is percent shade on the animal. 



5 Diffuse components (sky and substrate) of solar radiation are then calculated as 

follows: 

QSOBJ = ABSAN * Fabush * Area * PCTDIF * Qsolar,hor * (100. - ashade) / 
100 

QSSKY = ABSAN * Fa,sky * Area * PCTDIF * Qsolar ? hor * (100. - ashade) / 
10 100. 

QSRSB = ABSAN * Fa,grd * Area * (1.0 - ABSSB) * Qsolar,hor * (100. - 

ashade) / 100. 
QSDIFF = QSSKY + QSRSB + QSOBJ 

where, 

QSOBJ is incoming diffuse radiation reflected by an object nearby, 
Fabush is the configuration factor between the animal and the nearby 
object, 

Qsolar,hor is the incoming solar radiation on a horizontal surface, 
QSSKY is the diffuse solar radiation from the sky (i.e. solar radiation 

scattered by molecules in the atmosphere and by clouds), 
Fa,sky is the configuration factor between the animal and the 

skyQSRSB is the incoming solar radiation reflected from the 
substrate (i.e. ground or other surface), 
Fa,grd is the configuration factor between the animal and the ground, 
ABSSB is the absorptivity of the substrate, and 
QSDIFF is the total amount of diffuse solar radiation. 




Total solar energy absorbed (Qsolar,a) is thencalculated as the sum of energy 
30 absorbed on the top (i.e. dorsal; Qdorsl) and bottom (i.e. ventral; Qventr) parts of animal 510. 
Qdorsl = (QSDIR + QSOBJ + QSSKY) 
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Qventr = QSRSB 
Qsolar,a = Qdorsl + Qventr 



2. Estimating the metabolic heat generation needed to maintain the 
5 animal's core temperature (Qgen) 520 and computing the temperature of 

the skin (Tskin) 530. 

An estimate for the metabolic heat generation needed to maintain the animal's core 
temperature (Qgen) is made 520. Before calling the iterative guessing routine, Zbrent, the 
iterative guessing routine two bounding guesses are needed to start the process. This is done 

10 by first guessing for total metabolic rate (W) using the standard regression value in the 

literature for an animal at 28 C air and radiant temperature, humidity approximately 5%, and 
low wind speed, approximately 0.1 m/s. Humidity and wind speed were never reported for 
the experiments. MASS is in kg. The regression units are not J/s, so conversion is required. 
Qgen=Xl=(70.* MASS**0.75)*(4.185/(24.*3.6)) 

1 5 A second guess is created by arbitrarily adding 5 to the value of XI just calculated. A 
stability check is done on the value of XI to make sure it does not equal zero. Then a 
bracketing subroutine, ZBRAC, (Press, W. H, Flannery, B .P., Teukolsky, S.A., and 
Vetterling, W. T.. 1987. Numerical Recipes. The Art of Scientific Computing. Cambridge 
University Press, Cambridge ) is called. 

20 CALL ZBRAC(FUN,X1,X2,SUCCES) 

FUN is the name of the function containing the energy balance equation, XI and X2 are two 
starting guesses, SUCCES returns a value indicating success of the process. 
A criteria for success, TOL, in balancing the energy balance equation must be specified. 
Ideally TOL equals zero if heat fluxes in and generated in exactly equals heat fluxes out. 

25 In the program TOL = 0.01 *mass 

Then a solution is sought by invoking the energy balance equation via ZBRENT (Press et al., 
ibid.) and FUN. 

Qgen = ZBRENT(FUN,Xl,X2,TOL) 
Inside FUN the estimated value for Qgen is then equated to the value for the heat generation 

30 per unit volume (g ), 

g=Qgen/VOL 
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Respiratory water loss as a function of oxygen consumption is then done based on the heat 
generation requirement. No negative metabolic rates are allowed. 
If(Qgen.gt. 0.000) then 
CALL RESP 

5 else (Negative metabolic rate. No physiological meaning - set dummy values 

this guess) 
Qresp = 0.00000 
Qgen = -10.00000 
endif 

1 0 Net internal heat generation is calculated as 
Qgenet = Qgen - Qresp 
Net internal heat generation/unit volume is then calculated. 
Gn = Qgenet/Vol 

This is used for estimating skin temperature as in Eq. (53) (Porter et ai } 1994, at p. 162, 
1 5 herein incorporated by reference): 

Gn = g - [(n K) / R 2 ] * (Tcore - Tskin) 
where, 

n is the geometry constant (varies with shape), 
K is the thermal conductivity for tissue, 
20 R is the radius of the animal, 

Tcore is the core temperature of the animal, and 
Tskin is the skin temperature of the animal. 

Rearranging the equation, the Tskin is solved for 530. The value for Tskin is then 
25 used in subsequent calculations. 

3. Computing heat fluxes between the animal and its environment 540. 

The heat fluxes occurring between the animal and its environment are calculated by 
solving for the amount of heat exchange across the insulation of the animal using a porous 
insulation model to compute both the amount of heat entering the animal and the heat losses 
30 541 from the animal. An appendage model is also employed in the case of large birds and 
mammals to further incorporate heat fluxes across an animal's appendages. 
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a.) The porous insulation model. 

The porous insulation model of the present invention is run to calculate the amount of 
heat exchange across the insulation of the animal (Qin - Qout). For birds and mammals 
greater than 1 00 kg in size, an appendage model (see below) is also employed to more 
5 accurately calculate the heat exchange over appendage surface areas. 

Unlike the prior version of the porous insulation model [W.P. Porter et al 9 1994. 
Endotherm energetics; from a scalable individual-based model to ecological applications, 
Aust J. Zool. 42:125-62, Appendix equations (24) to (53) on pp. 159-162 and description of 
the model and methods on pp. 127- 137, the disclosure of which is incorporated herein by 

10 reference], the porous insulation model of the present invention incorporates a solar radiation 
calculation to broaden its applicability to daylight conditions. In this way, the porous 
insulation model enables a user to calculate the amount of Qin including both when the 
animal is exposed to sunlight and when the animal is in the dark. 

Solar radiation (Qsolar) is incorporated into the porous insulation model by assuming 

1 5 that solar radiation is absorbed very close to the fur/feather-air interface, which is usually the 
case for bird feathers and dark, dense fur. Absorbed solar radiation (Qsolar,a) heats the fiber 
elements, which then emit infrared radiation outward (QIR,out) toward the sky and inward 
(QIR,in) through the porous insulation. The watts of absorbed solar radiation (Qsolar,a) are 
treated as an additional source of thermal radiation from the sky for the half of the animal 

20 exposed to the sky. Thus, the diffuse IR equations already in the prior model were also used 
for incorporating absorbed solar radiation in the model of the present invention. To do this, 
an additional term for absorbed solar radiation was added to equation (52) of the original 
model (Porter et al 9 1994, eq. (52) at p. 162). 

Computing IR fluxes in. The Swinbank formula is applied to solve for the clear sky 

25 temperature (Tsky) in degrees Kelvin, based on the air temperature at 2 m in height (Tair,2m) 
in degrees centigrade. 

Tsky = 0.0552 * (Tair,2m + 273.)**(1. 5-273.) 

30 IR fluxes from other objects overhead of the animal (i.e. vegetation) that are 

'substituting' for sky, are also calculated. The temperature of vegetation (Tveg) is set to the 
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Tair,2m to solve for the IR fluxes from sky (QIR,sky } and vegetation (QIR,veg } : 
Tveg =Tair,2m 

QIR,sky = C3 * ( Tsky + 273.1 5)* *4 
QIR,veg = C6 * ( Tveg +273.1 5)* *4 

Next, the porous media model checks for whether fur or feathers are optically thick or 
not. If optically transparent [i.e. no effective fur or feathers or no insulation present (as is the 
case for many ectotherms)], it is assumed that no solar radiation is absorbed by the fur or 
feathers. If, on the other hand, the fur or feathers are optically thick, solar radiation is treated 
as though it is absorbed close to the fur-air interface and is treated as an IR source. 

IF(B1ARA(1) .It. 0.001), then fur/feathers optically transparent (no effective fur) 

and, 

QIR,above = QIR,sky + QIR,veg 
QIR,below = sig * ( Tgrd(i) + 273.1 5)* *4 
where, 

BlARA(l) is a variable that indicates the optical thickness of the fur 
or feathers; it is the exponent of Beer's Law {i.e. the extinction 
coefficient (x)times the depth of fur/feathers (t)), 
QIR,above is the IR fluxes from above the animal, 
QIR,below is the IR fluxes from below the animal, 
sig is the Stefan-Boltzmann constant, and 
Tgrd(i) is the ground temperature at location (i). 

Else, optically thick, fur present and solar radiation absorbed, qdorsl from 

above and qventr (solar reflected from the ground) is treated as diffuse 

radiation absorbed by the fur, 

QIR,above = QIR,sky + QIR,veg + qdorsl 

Tsky = (QIR,sky / (C3+C6))**0.25 - 273.15 

QIR,below = sig * ( Tgrd + 273.1 5)* *4 + qventr 
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Tgrd = (QIR,grd / sig)**0.25 - 273.15 



Whichever the case, the amount of radiant heat entering the animal from external 
sources (QIR,in) is: 

QIR,in = QIR,above + QIR,below 
as illustrated_ in FIG.3a. 

Computing heat losses via conduction, convection, evaporation and IR flux. 541. 
Equation (52) of the original model (Porter et al, 1994. eq. (52) at p. 162, see PP. 156-169 
for explanation of variables) is: 

CI = - [4a ((Tskin + Tsky) / 2) 3 (Tskin - Tsky)(l - 2E 3 (^ T Z L )) + h c (Tskin - Tair) + 

2c(T 4 skin - tsky)Ei(fi T Z L )] I [(k eff + k r ) + (h c + 2o((Tskin + Tsky) / 2) 3 (1-2 
E 3 (J3 T Z L )))Z L ] 

In the new version (Porter et. al, 2000), the equation has been modified so that the 
sky, the ground, and nearby objects (bush, boulder, etc.) are treated separately and then 
combined. The equations have been normalized for more compact form. 

The sky portion of radiant heat exchange (TAVsky) is, 

TAVsky = ((Tsky + Tskin) / 2.) + 273.15 

KRADSKY = (16.0 * sig * TAVsky**3.) / (3. * BETARA(2)) 

HR = Fa,sky * 4. * EMIS * sig * (Tsky + 273.15)**3 

BIR = HR * ZFURAR(2) / (KEFARA(2) + KRADSKY) 

BICV = HC * ZFURAR(2) / (KEFARA(2) + KRADSKY) 

Dl = 1. + BICV + 2. * BIR * (0.5 - (0.333 - E4B1 / B1ARA(2))) 

C1N = (BICV * (1. - Tkair / TKskin) + BIR * (1 . - Tksky / TKskin)) 

Cldors = ClN/Dl 

PSI = 0.854 + (ALOG10 (BICV) * 0.149) 

QRADSK = PSI * Areaskin * Fa,sky * KRADSKY * (CI dors * TKskin / 
ZFURAR(2)) 
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where, 

TAVsky is the average temperature on the dorsal half of the animal 
(degrees K), 

KRADSKY is the fur/feather radiant conductivity in the sky direction, 
5 BETARA(2) is the IR extinction coefficient of fur/feathers, 

HR is the radiant heat transfer coefficient, 

EMIS is the emissivity of the animal, 

BIR is a dimensionless variable for IR heat transfer, 

BICV is a dimensionless variable for convection heat transfer, 
1 0 ZFURAR is the depth of fur, 

KEFAR is the effective thermal conductivity of fibers and air for 
conduction, 

HC is the convection heat transfer coefficient, 

Dl is dimensionless, 
15 C1N (CI normalized) is dimensionless, 

PSI is a correction factor, 

TKair,TKskin, TKsky are the temperatures of air, skin, sky in degrees 
K, 

QRADSK is the incoming solar IR radiation from the sky through the 
20 porous media, 

Tskinis the skin temperature of the animal, and 

Subscript notation: (1) = overall, (2) = dorsal, (3) = ventral fur 

The following calculations continue in the same way as those above, but for a bush or 
25 nearby large object portion of radiant heat exchange, instead of for sky (definitions above are 
easily substituted for variables below with respect to bush instead of sky and will not be 
defined separately): 

TAVbsh - ((Tbush + Tskin) / 2.) + 273.15 
30 KRADBUSH = (16.0 * SBK * TAVbsh **3.) / (3. * BETARA(2)) 

HR - Fabush * 4.* EMIS * SBK * (Tbush + 273.15)**3 
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BIR = HR * ZFURAR(2) / (KEFARA(2) + KRADBUSH) 
BICV = HC * ZFURAR(2) / (KEFARA(2) + KRADBUSH) 
Dl = 1. + BICV + 2.* BIR * (0.5 - (0.333 - E4B1 / B1ARA(2))) 
C1N = (BICV * (1 . - TKair / TKskin) + BIR * (1 . - Tkbush / TKskin)) 
Clbush = ClN/Dl 

PSI = 0.854 + (ALOG10 (BICV) * 0.149) 

QRADBU = PSI * Areaskin * Fabush * KRADBUSH * (Clbush * Tkskin / 
ZFURAR(2)) 

The ground portion of radiant heat exchange (ventral part of animal), 

TAVgrd = ((TGRD + Tskin) / 2.) + 273 . 1 5 
ktop = 16.0 * SBK * TAVgrd**3. 
kbottom = 3. * BETARA(3) 
KRADGRD = ktop / kbottom 

HR = Fa,grd * 4.* EMIS * SBK * (((T2 + Tskin) / 2.) + 273.15)**3 

BIR = HR * ZFURAR(3) / (KEFARA(3) + KRAD) 

BICV = HC * ZFURAR(3) / (KEFARA(3) + KRAD) 

Dl = 1.+ BICV + 2.* BIR * (0.5 - (0.333 - E4B1 / B1ARA(3))) 

C IN = (BICV * (1 . - TKair / TKskin) + BIR * (1 . - TKgrd / TKskin)) 

Clvent = ClN/Dl 

PSI = 0.854 + (ALOG10 (BICV) * 0.149) 

QRADGR = PSI * Areaskin * Fa,grd * KRADGRD * (Clvent * Tkskin / 
ZFURAR(3)) 

The net radiant heat exchange (QRAD) is then the net radiant fluxes between the 
animal and the sky, ground and bushes (or other objects), i.e. QRAD = QIR,in - QIR,out (in 
FIG3a). 

QRAD = QRADSK + QRADGR + QRADBU 
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Conduction not lost by outward radiation from hairs within the fur equals convection 
at fur-air interface (Qconv). 

Dorsal conduction (QcondD) (via air (Qcond,air) and fur/feathers(Qcond,hair)): 

QcondD - PSI * (Areaskin / 2.) * KEFARA(2) * (CI dors * TKskin / 
ZFURAR(2)) 

Ventral conduction (QcondV) (via air (Qcond,air) and fur/feathers(Qcond,hair)): 

QcondV = PSI * (Areaskin / 2.) * KEFARA(3) * (Clvent * TKskin / 
ZFURAR(3)) 

Total heat loss by conduction (Qcond) is then: 

Qcond = QcondD + QcondV 

CI = (Cldors + Clvent)/2. 
The fur energy flux (Qfur) is the total heat loss from the skin and fur. 

Qfur = Qcond + QIR ? out + Qevap,skin + Qconv 
The whole animal energy balance at the skin is then, 

Qnet = Qin + Qgen - Qout 

where, 

Qnet is the net heat exchange of the whole animal; 
Qin is the total energy entering the animal = Qsolar,a + QIR,in; 
Qgen is the metabolic heat generation needed to maintain the animal's 
core temperature (estimated earlier in the process); and, 
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Qout is the total energy leaving the animal = Qfur + Qevap,resp (where 
Qevap,resp is calculated below by the lung model as part of the 
mass balance analysis). 

5 Computing the ERROR by solving the heat balance equation 550. Assuming that the 

value of Qnet will approach 0 when the value of Qgen is accurate, the ERROR is substituted 
for Qnet in the equation and solved for as follows, 

ERROR = Qin + Qgen - Qout 

10 

The actual value of ERROR is then compared to the user-specified maximum 
acceptable value for Error, TOL, 555. If the actual ERROR is greater than the specified 
value, a new estimate for Qgen is generated 520 and a new value for Tskin is calculated 530 
for subsequent use in the following computations of heat losses 540. 

1 5 If the ERROR is less than or equal to the specified maximum value, the value of Qgen 

is deemed to be accurate, and the value for Tskin previously calculated on its basis 530 is also 
deemed accurate. In that case, the thermoregulatory model 537 is actuated (see below) and a 
calculation of the core-skin temperature gradient (Tcore-Tskin) is made and compared to a 
user-specified value (for example, 0.5 degrees C). 

20 If the core-skin gradient is less than 0.5 C, the animal must respond behaviorally or 

physiologically to enable proper cooling. In this case, the thermoregulatory model 537 
invokes the behavioral and physiological options in a predetermined order. After an option is 
chosen, the animal model 500 proceeds to estimate the Qgen 520 and proceed through the 
aforementioned steps until evaluating the core-skin temperature gradient again. If the 

25 gradient is still not satisfied, the next option is set and calculations repeated until the 
minimum core-skin gradient is satisfied. 

When the core-skin temperature gradient is satisfied, (i.e., it is greater than or equal to 
0.5 C), the animal is assumed to be able to cool itself properly. The animal model 500 then 
proceeds to begin another round of calculation over every other hour of a day to calculate an 

30 average day of each month, then over every month to eventually calculate annual values for 
the various energetic variables. 
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b. ) The appendage model. 

An appendage model was also added to the animal model 500 to more accurately 
calculate the amount of heat entering and lost from the animal via its appendages when the 
model is applied to large birds (such as, for example, the rhea, the cassowary, and the ostrich) 
and to mammals of a size greater than or equal to 1 00 kg. As applied to those sorts of 
animals, the appendages are large and may constitute a significant percentage of the surface 
area of the standing animal over which heat may be gained or lost. The geometries used to 
model a large standing bird (cylinders and ellipsoid) compared to a sitting bird (ellipsoid) 
vary accordingly. But without porous insulation and respiratory water loss, the appendage 
model is used to calculate heat loss in the radial dimension from the appendages of an animal. 

c. ) The thermoregulatory model 537. 

The thermoregulatory model yields modified microclimate/physiology conditions that 
allow for a stable solution for body temperature when there is heat stress (i.e., when the core- 
skin temperature gradient (Tcore-Tskin) is less than 0.5 degrees C, as discussed above). The 
thermoregulatory model 537 solves for the behavioral and/or physiological changes needed to 
maintain the animal's core temperature in heat stress conditions. 

Options in a behavioral and physiological options hierarchy are chosen sequentially. 
Initially, the animal model 500 begins by assuming the animal is in a "null" condition. The 
"null" condition is user-specified and represents a circumstance where the animal is exposed 
to the sun and sky as much as possible given user specified constraints (e.g., regarding 
whether the animal is diurnal, nocturnal or crepuscular; whether it is arboreal or fossorial; the 
percent shade available and degree of cloud cover in the animal's environment; and the like). 

If the core/skin gradient is not greater than or equal to the specified minimum (e.g., 0.5 
degrees C), then the thermoregulatory model's 537 options hierarchy is invoked. In the first 
level, the animal may proceed to modify its behavior by, for example, seeking shade on the 
ground, climbing to cooler temperatures in the shade, retreating below ground, or acting 
otherwise depending on the user specified behavioral options. Option-dependent values for 
variables are set and the animal model 500 proceeds to begin again to estimate Qgen 520 and 
following (as described above, see FIG. 2) until the core-skin gradient is once again 
evaluated. If the chosen option is insufficient to maintain at least the minimum core-skin 
gradient, then the next behavioral option is exercised, and so on through the behavioral 
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options until the minimum core/skin temperature gradient is satisfied. 

If the behavioral options fail to satisfy the minimum core-skin temperature gradient, 
then the animal is allowed to exercise physiological options and the model iterations repeated 
as with the behavioral options. The user may specify that the animal first allow its core 
temperature to rise before invoking evaporative heat loss, or alternatively, the user can specify 
that the animal invoke evaporative heat loss first and, failing that, allow its core temperature 
to rise. A core temperature increase will always allow for a stable solution and may be 
exercised as a last option. 

Many different and varied combinations of behavioral and physiological options are 
possible and specifiable by the user. Each will generate its own set of energetic costs that will 
subsequently feed into the animal model's 500 calculation of the amount of discretionary 
energy and water. The thermoregulatory model's 537 evaluation of the behavioral and 
physiological options available to an animal may yield significant insights into how best to 
manage the environment of an animal for minimal stress, optimal activity time and other 
related effects. 

ii. Performing the mass balance analysis 502 

The mass balance analysis of the present invention solves for the amount of 
discretionary mass from food (Mf,discr) [given the amount of food needed to meet the 
metabolic heat generation required to maintain the animal's core temperature (Qgen)], the 
amount of discretionary mass from water (Mw,discr), and the amount of discretionary energy 
(Qdiscr) available to the animal 

The mass balance equations of the animal model 500 generally solve for the total 
discretionary mass from food and water (Mdiscr) by calculating the mass of food and water 
absorbed by the animal (Mf,a, Mw,a) [equaling the mass of food and water intake by the 
animal (Mf,i, Mw,i) minus the mass of food and water loss by excretion and evaporation 
(Mf,d, Mw,u, Mw,evap)], and the mass of food and water metabolized by the animal (Mf,met, 
Mw, met). To do this, the animal model 500 employs molar balance models for the lungs 570 
and gut 560 to allow for the coupling of food ingested and/or absorbed to respiratory gas 
exchanges and evaporative water loss from the respiratory system. 

One of the improvements in the animal model 500 is the addition of the distributed 
respiratory water loss, which represents lungs that span most of the body cavity. This 
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innovation gives much better agreement of predicted metabolic rates with measured values. 
A skin subroutine is also used to calculate evaporative water loss from skin (Mw,evap,skin). 

1. Solving for the amount of discretionary mass from food (Mf,discr). 
5 a.) The gut model 560. 

The gut model 560 couples food ingested to respiratory gas exchanges and evaporative 
water loss from the respiratory system. FIG. 4 shows the system diagram for the molar 
balance gut model 560. It is related to the well-known batch reactor and plug flow model 
originally developed in chemical engineering and subsequently applied to animal digestive 
10 systems (see, Penry, D.L. and P. A. Jumars, 1987. Modeling animal guts as chemical 

reactors. Am. Nat. 129:69-96). The model used in the present invention allows for any type 
of ingested food consisting of user specified proportions of carbohydrates, lipids, proteins and 
water content. The food can enter the gut any time during activity time in any amount, 
U subject to the constraint that the volume of food ingested per day may not exceed the wet 
yp 15 mass of the animal. The energy value of absorbed carbohydrates, lipids, and proteins is well 

known {see, Schmidt-Nielsen, K. 1979. Animal physiology, adaptation and environment. 2d 
M ed. ? Cambridge Univ. Press, Cambridge, MA). 

fn The gut model 560 determines an upper bound for absorbed food. It does not 

!l explicitly consider important constraints on food intake that are affected by how food is 
£020 distributed in the environment, which affects encounter probability, and what the handling 
5y time is, which is a function of feeding apparatus morphology. Both encounter probability and 
p handling time may affect how much food is actually ingested. The present version of the gut 
model waits until a 24 hour activity simulation cycle starting at midnight has been completed. 
Then it looks back at the day's energy requirements and uses that to calculate how much of 
25 the food would have to be processed in the gut that day to meet those demands. Other cycle 
lengths could be used instead, depending on the user's interests. 

The user specifies the food properties. They include percent of proteins, 
carbohydrates, lipid, and percent dry mass relative to wet mass. There is an upper bound 
constraint that the animal may never consume more wet weight food per day than its body 
30 mass. 
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The gut model 560 deals only with the part of the gut associated with the actual 
digestive process. Handling time and internal storage of food in the stomach(s) is not 
considered. These will place additional constraints on gut function. Because retention time 
may vary from approximately one-third of a day to more than a day, depending on body size, 
a day's energy requirements and the mass flow of food per day needed to meet those 
requirements are used. The total energy requirements are set to some factor of the energy 
needed to maintain core temperature on any particular day simulated. For example, the total 
energy requirements may be set to equal 3.5 times the energy needed to maintain the animal's 
core temperature on a particular day. 

A digestive efficiency is set to a value depending on food type. For example, a 
digestive efficiency of 57% has been measured for seeds in small rodents and is very close to 
the value for digestive efficiency by cattle (61%). The model requires empirical data on 
digestive efficiency and food properties as defined above (these data are typically provided by 
the user but may instead be provided by the animal characteristics database 190 if appropriate 
data are contained therein). The gut model 560 assumes that the gut can alter flow rates to 
maintain the same digestive efficiency for a particular type of food. 

Animals typically do not overeat or gain weight unnecessarily. The gut model 560 
assumes that the amount of food ingested is determined by today's metabolic needs and is 
constrained by maximum volumetric flow rate less than or equal to body mass. The user 
defines percent digestive efficiency, percent water in feces, percent urea in urine, and a 
multiplier from 1 to 7 (timbas) above basal metabolic energy requirements to specify food 
intake to meet current needs. The user also defines percent protein, percent fat, percent 
carbohydrate, and percent dry matter in ingested food. 

The dry grams of food needed to meet the metabolic demand are estimated based on 
food properties and the constraints specified above. 

The gut model 560 is related to work by Penry and Jumars (1987, ref. above, 
incorporated herein by reference). Their equation (9), p. 73, for a batch reactor model is: 
time = C A 0 Jo Xaf (d XA / -r A) 
where 

XA = mole fraction of species A (dimensionless) 
XA = (CA 0 - CA) / CA 0 
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CA = food component, A, concentration (moles) 

CA 0 = initial concentration of food component, A ? in moles. 

Many of the components of reactor models are not known. Since the endotherm 
5 thermoregulatory model (described in previous section) can calculate metabolic heat 

production to maintain a specified core temperature (Qgen), it is straightforward to calculate 
the grams dry weight, and then wet weight of user-defined food to meet the metabolic 
demand. The following pseudo-equations are an example of these calculations for a seed diet 
by way of illustrating the method. 
1 0 The amount of fat, protein, and carbohydrate per gram dry mass of food (gfatpg, 

gprotpg, gcarbpg, respectively) is first computed from user specified percentages of these 
components of food (pctfat, pctpro, pctcar, respectively), as in: 

gfatpg = pctfat (g fat / g dry mass) 
1 5 gprotpg = pctpro (g protein / g dry mass) 

gcarbpg = pctcar (g digestible carbohydrate / g dry mass) 
gsum = gfatpg + gprotpg + gcarbpg 

Next undigested mass per gram dry mass of food (gundig) is calculated and used to 
20 compute total carbohydrates (TOTCARB): 

gundig = 1.00 - gsum 
TOTCARB = gcarbpg + gundig 

25 Then joules per gram dry food (TotJpgram) is calculated: 

fatJpg = gfatpg * (9400. * 4.185) (g fat / g dry food * calories / g fat * 

J/calorie) 
protJpg = gprotpg * (4199 * 4.185) 
30 carbJpg = TOTCARB * (4200. * 4.185) 

TotJpgram = fat Jpg + protJpg + carbJpg 
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Not implemented here is the variable protein digestion efficiency (DigEff ) for 
ruminants versus monogasters (approximately 75% vs. 13%). While proteins are important 
from a protein and amino acid balance standpoint, in terms of absolute energy available, it is a 
5 small factor, since alfalfa (reference number 100056) vs. grass forage (reference number 
102250) varies from approximately 8.6% protein to 18.1% protein (US-Canadian nutrient 
compound composition tables, 1994). Seed protein is 15% (see, Kelrick, M.R. and J. A. 
MacMahon, 1985. Nutritional and physical attributes of seeds of some common sagebrush- 
steppe plants: some implications for ecological theory and management. J. Range Manage. 
10 38(l):65-69). 

Calculation of joules absorbed per gram of dry food (Jabspgr) is: 



Jabspgr = fatJpg + protJpg + DigEfPcarbJpg 



1 5 The moles available/g dry food (Totmolpgram) can now be calculated: 

Totmolpgram = gfatpg * 850 + gprotpg * 137 + gcarbpg * 180 

The joules of chemical energy needed to be absorbed by the gut /day (Food J) = joules 
20 per second (J/s) * seconds/hour (s/h) * 24 hours/day * basal multiplier for metabolic rate (i.e. 9 
activity above resting; timbas), can then be calculated: 

FoodJ = Qmetab J/s * 3600 s/h * 24 h/d * timbas 

25 The dry mass needed to be absorbed per day to meet the metabolic demand is 

then, 

DRYA - FoodJ/Jabspgr 



The grams of dry food that must be ingested per day (Drymas) and then the total 
30 moles of each of the foodstuffs (Totmoles) are calculated: 
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Drymas = DRYA / (1 .0 - gundigest) 

Total moles available for absorption = moles/dry gram * dry grams 
Totmoles = (gfatpg * 850. + gprotpg *137 + gcarbpg *180.) * drymas 

The joules of food needed in the gut per day to meet resting metabolism plus user 
specified activity level = (J/DAY)*basal multiplier: 

FoodJ = DaysMETAB * timbas 
The needed grams of food absorbed per day = J needed/day / Food J/g : 

DRYABS = FoodJ / Jabspgr 

The grams dry food ingested per day = g absorb/DigEff: 

Drymas = DRYABS / (1.0 - gundigest) 
Wetmas = Drymas / PctDry 

The program band tests for whether the required food mass per day exceeds body 
mass. If it does, absorbed mass is reset to the maximum value allowed by body weight and 
energy available is adjusted accordingly. 

Upprlim = Gmass 

If (Wetmas .greaterthan. Upprlim), then 
Wetmas = Upprlim 
Drymas = Wetmas * PctDry 
DRYABS = Drymas * DigEff 
Jpdaymax = DRYABS / Jabspgr 
Jpsavail = Jpdaymax / (24.* 3600.) 

Endif 
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Once it is known that the upper bound has not been exceeded, the day's joules of 
energy absorbed (DaysJabs) can be calculated. 

DaysJabs - Jabspgr * DRYABS 

The energy available, Avlnrg, in joules/s available can be calculated 
Avlnrg = DaysJabs / (24.*3600.) 

The g wet weight through the gut/day is then calculated (CRTOT) to monitor part of 
the water balance of the animal: 

GW = (Drymas - PCTDRY * Drymas) / PCTDRY 
CRTOT = gw + Drymas 

The hourly dry grams absorbed of protein, fat and carbohydrate (GPROTH, GFATH, 
GCARBH, respectively) are also tracked for use in the growth model (see, Adolph and Porter, 
1996. ref. above). 

Dryinphr - DRYABS / 24. 
GPROTH = Dryinphr * gprotpg 
GFATH = Dryinphr * gfatpg 
GCARBH = Dryinphr * gcarbpg 

Size dependent constraints appear when the energy (mass) needs for the day exceed 
the daily gut volume flow rate for the body size. Available activity time available is 
calculated. That can be used to back calculate required handling time of the food needed to 
meet the energetic demands and to estimate time efficiency of mass/energy acquisition for all 
body sizes for a given food type and configuration. 

The model assumes 0.40 g water/g protein oxidized (from, Hainsworth, F.R. 1981. 
Animal physiology. Addison- Wesley Publ. Co., Reading, MA). 
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Proteins, The approximate gram molecular weight of amino acids is 137 g/mole 
based on an average over all 22 amino acids (based on on information from the CRC 
Handbook of Chemistry & Physics (1996)). 

Lipids, Data from Guyton (Guyton, A.C., 1991 . Textbook of medical physiology. 8 th 
5 ed., W.B. Saunders, Philadelphia, PA) are used. The gut model assumes that triglycerides are 
used for energy. They are stearic acid (880 g/mol), oleic acid (879 g/mol) and palmitic acid 
(754 g/mol). Based on these data the model assumes as an average 850 g/mol. The model 
uses 4.7 calories/ml 0 2 which is 9400 calories/g fat (see, Kleiber, 1975. The fire of life. An 
introduction to animal energetics, 2d ed. Krieger Publ. Co., Huntington, NY). There are 1.07 
1 0 g water produced/g lipid oxidized (see, Hainsworth 1981, ref. above). 

Carbohydrates, The model assumes that glucose data values are 180 g/mol; 5.0 cal/ml 
0 2 ; 4200 calories/g (see, Kleiber 1975, ref. above) and that there are 0.56 g water produced/g 
carbohydrate oxidized (see, Hainsworth 1981, ref. above). 
% Digestion capability data used by the model come from several sources (see, Weiner, 

S 15 J. 1 987. Limits to energy budget and tactics in energy investments during reproduction in the 
Djungarian hamster, Phodopu sungorus sungorus, Pallas 1770, in Reproductive energetics in 
^ mammals. Symp. Zool. Soc. London 57:167-187; Kenagy, G.J. et al. 1989. Energy 
ffi requirements for lactationa nd postnatal growth in captive golden mantled ground squirrels. 
q Physiol. Zool. 62:470-487; Peterson, C.C. et al. 1990. Sustained metabolic scope. Proc. 
ffl20 Natl. Acad. Sci. 87(6):2324-2328). Seed component values come from Kelrick and 
hi MacMahon (1985, ref above). 

£7 2. Solving for the amount of discretionary mass from water (Mw,discr). 

Calculating the discretionary mass from water (Mw,discr) requires running the gut 
model 560 to calculate the mass of water absorbed (Mw,a). The lung model 570 calculates 

25 respiratory water loss and the skin evaporation model (i.e. the skin subroutine) computes 
water lost by evaporation due to convection. The difference between water absorbed and 
water loss is the discretionary water (Mw,discr). Water excreted in urine (Mw ? u) is 
calculated based on the amount of protein in the diet (see below). As mentioned above, an 
improvement in the animal model 500 of the present invention is the addition of the 

30 distributed respiratory water loss, which represents lungs that span most of the body cavity. 
This innovation gives much better agreement of predicted metabolic rates with measured 
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values and is achieved by use of the lung model 570. 

a.) The lung model 570. 

This subroutine is a molar balance for computing water loss from breathing. It uses 
the oxygen demand for maintaining a user-specified core temperature to compute the amount 
5 of air in and out of the lungs. The temperature of the exhaled air relative to ambient air 
temperature is assigned by the user. This is typically anywhere from 1-3 ° C higher than 
ambient air. This in effect allows for recovery of water used to saturate air in the lungs, 
because air leaving the nares is assumed to be saturated at the exhaled air temperature. 
Below are the basic calculations performed by the lung model 570. 
1 0 The partial pressure of oxygen at current conditions is, 

P0 2 = BarPrs * PCT0 2 
where, 

BarPrs is barometric pressure (pascals) 
1 5 PCT0 2 is the percent 0 2 in air 

P0 2 is the partial pressure of 0 2 (pascals) 

The partial pressure of oxygen at standard temperature and pressure (STP) is 

20 RefP0 2 = 101325 * Rpct0 2 

The lung model 570 calculates the total heat production needed to maintain a user- 
specified core temperature (Qgen), given the current environmental conditions and the 
animal's characteristics. Oxygen consumption is based on the calculated heat generation 
25 estimate to maintain body temperature (Qgen), corrected for substrate utilized and is 
calculated as follows: 

Liters of 0 2 per second at STP: (using data for equivalencies from Kleiber 
1975, see ref. above) 



30 



if(RQ .eq. 1.0) then, 
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carbohydrates are worth 4200 cal/g 

liter(STP)/s = J/s / ((cal/J) * (kcal/cal) * (liters 0 2 /kcal)) 

0 2 STP = Qgen * Timbas / (4.185 * 1000 * 5.0) 

else 

if(RQ .le. 0.7) then, 

fat; fat's worth 9400 cal/g 

0 2 STP = Qgen * Timbas / (4.185 * 1000 * 4.7) 

else 

protein (RQ=0.8); proteins worth 4300 cal/g on average 
0 2 STP - Qgen * Timbas / (4.185 * 1000.* 4.5) 

where, 

Timbas is the basal multiplier 

Qgen is the current guess for total heat generation 

RQ is the respiratory quotient 

Converting O2STP to the volume of O2 at the animal's core temperature, the general 
equation (V1*P1/T1) = (V2 * P2 / T2) is rearranged to solve for volume of oxygen 
consumed and the ideal gas law is used: 

n = PV/RT (Ideal gas law: number of moles from pressure, volume, 
temperature) 

V0 2 CON = (0 2 STP * RefP0 2 / 273.15) * ((Tlung + 273.15) / P0 2 ) 
where, 

O2MOLC - 0 2 consumed / second - P0 2 * V0 2 CON / (RGC * (Tlung 
+ 273.15)) 

Tlung = the temperature in the lungs obtained from a solution for the 
average temperature from the core to the skin obtained from 
integration. 

MOLES/s O2, N2, & dry air AT 1 : (1 refers to entrance of respiratory system, 
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while 2 refers to the exit, see FIG 4) (Air flow = f(0 2 consumption) 
0 2 MOLl = O2MOLC / (EXTREF / 100) 
N 2 MOLl = 0 2 MOLl * (PCTN 2 / PCTO2) 
Demand for air = f (%02 in the air and elevation) 

Note that as long as all 3 percentages add to 100%, no change in air flow occurs, 
unless one corrects for change in %0 2 in the air and elevation changes relative to sea level. 

AIRML1 = 0 2 MOLl * ((PCTN 2 + PCT0 2 + PCTC0 2 ) / PCT0 2 ) * (Rpct0 2 / Pct0 2 ) 

(RefP0 2 /P0 2 ) 
Air volume at STP (liters/s) * (s/hour) 
AirVol = (Airmll* RGC * 273.15 / 101325.) * 3600. 

Computing the vapor pressure at saturation for the subsequent calculation of actual 
moles of water based on actual relative humidity, 

rhsat= 100 

CALL WETair(Tair,WB,rhsat,dP,BarPrs,E,ESAT,VD, 
RW,TVIR,TVINC,dENAIR,CP,WTRPOT) 

Moles water/s in at 1 (entrance, see FIG. 4) based on relative humidity. Note that 
humidity is set to 99% in Main if animal is in burrow for the current hour. 

WMOL1 = AIRML1 * (ESAT * (RELHUM / 100.)) / (BarPrs - ESAT * 

(RELHUM/ 100)) 
MOLES/s of dry air AT 2: (EXIT) 
0 2 MOL2 = 0 2 MOLl - 0 2 MOLC 
N 2 MOL2 = N 2 MOLl 
C0 2 MOL = RQ * O2MOLC 

Total moles of air at 2 (exit, see FIG. 4) will be approximately the same as at 1 
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(entrance), since the moles of O2 removed = approx. the number of moles of CO2 added. 
Avogadro's number specifies the number of molecules/mole. 



AIRML2 = O2MOL2 + CO2MOL) * ((PCTN 2 + PCTO2) / PCT0 2 ) * (Rpct0 2 / 
Pct0 2 ) * (RefP0 2 /P0 2 ) 

Setting up call to WETair; TEMP. OF EXHALED AIR close to AIR TEMP, assuming 
saturated air at exhalation, 

RELXIT = 100 

Kangaroo physiology (Schmidt-Nielsen, K., et al. 1970. Counter-current heat 
exchange in the respiratory passages: Effect on water and heat balance. Resp. Physiol. 9: 
263-276) shows Ta,exit - Ta = 3C. 

TAEXIT = Tair + DELTAR 

Not letting internal air go below 0 C for stability purposes: 

If (Taexit.lt. 1.0) then 
Taexit= 1.0 

Endif 

Calculating saturation vapor pressure (Esat) at exit temperature. 

CALL WETair (TAEXIT, WB,repELXIT,dP,BarPrs,E, 

ESAT,VD,repW,TVIR,TVINC,dENAIR,CP,WTRPOT) 
WMOL2 = AIRML2*(ESAT/(BarPrs-ESAT)) 

ENTHALPY = U2-U1, internal energy only, i.e., LAT. HEAT OF VAP are 
only involved, since P,V, and T are assumed to be constant, thus not 
significant. 
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FLOW ENERGY, PV. (H = U + PV) 



Moles/s lost by breathing: EVPMOL = WMOL2-WMOL1 

5 Grams/s lost by breathing = moles lost * gram molecular weight of water: 

Mw,evap,resp = EVPMOL * 18. 

A cap is put on water loss for small animals in very cold conditions by assuming they 
1 0 will seek more moderate conditions if they exceed this cap. This will improve stability for the 
solution method (based on data from W.R. Welch. 1980. Evaporative Water Loss from 
Endotherms in Thermally and Hygrically Complex Environments: An Empirical Approach for 
Interspecific Comparisons. J. Comp. Physiol. 139: 135-143). The maximum value recorded 
jj for Prairie Dogs was 0.6 g/(kg-h) = 1 .667 x 1 0**-4 g / (kg-s). The highest recorded rate was 
jO 15 for a resting deer mouse at 8 g/kg-h = 2.22E-03 * (Edwards & Haines. 1 978. J. Comp. Physiol. 
§ 128: 177-184 in Welch, 1980). For a 0.01 kg animal, the max. rate would be 1.67 x 10**-6 



H20 



Tesval = 2.22E-03 * timbas * gmass / 1000 

If (gevap .gt. Tesval) then, 
gevap = Tesval 

endif 



25 Kg/s lost by breathing: 

gperhr = gevap * 3600 
KGEVAP = GEVAP / 1000 



30 



Latent heat of vaporization from sub. DRY AIR (J/kg) 
HTOVPR = 2.5012E+06 - 2.3787E+03 * Tlung 
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Heat loss by breathing (J/s) = (J/kg) * (kg/s) 
Qresp = HTOVPR * KGEVAP 

5 

b.) The skin subroutine. 

The equations used to solve for the mass of water lost by evaporation from skin 
(Wevap,skin), is 

Wevap,skin = Aeff * h D *(VD S - VD A ) 
10 Where Aeff is the effective surface area that behaves as a free water surface 

ho is the mass transfer coefficient computed from the heat transfer coefficient, 
h c , using the Colborn Analogy (Bird, R.B., Stewart, W.E., and E.N. Lightfoot 1960. 
Transport Phenomena. Wiley and Sons, New York) 

VDs is the water vapor density at saturation at skin temperature 
1 5 VDa is the water vapor density at saturation at air temperature outside the 

animal's boundary layer. 
The grams of free water absorbed from food (Mw,f) per hour, equals the total wet 
weight absorbed per hour minus the dry weight absorbed per hour, as in: 
Mw,f = (GTOTHR - PCTDRY * GTOTHR) / PCTDRY 

20 

The grams of metabolic water (Mw,i) formed per hour is, 

Mw,i = GPROTH * 0.40 + GFATH * 1.07 + GCARBH * 0.56 

The grams of water lost in urine produced (Mw,u) per hour is, 
25 GUREA = GPROTH * 0.343 

GURINE - GUREA / PTUREA 
Mw,u = (GURINE - GUREA) / 1 .0474 

The water in feces (Mw,d), is 
30 g wet weight absorbed/hr = DigEff * g processed/hr 

gproch = (GFREWA + GTOTHR) / DigEff 
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WETFEC = gproch - (GFREWA + GTOTHR) 
Mw,d - WETFEC * PFEWAT 



The total amount of discretionary water available (Mw,avail) per hour from gut and 
5 excretory system is calculated as a function of diet, 

Mw,avail = Mw,f + Mw,i - Mw,u - Mw,d 

Total discretionary mass from water is then, 

10 

Mw,discr = Mw,avail - Mw,evap,resp - Mw,evap,skin 

3. Solving for the total amount of discretionary energy (Qdiscr). 

The total amount of discretionary energy available to the animal (Qdiscr) is 
Ci 1 5 then calculated as follows: 

Qdiscr = Qavail - Qgen 

fjl iii. Calculating total activity time 

P~ The total annual activity time of the animal is a key variable linking the individual 

03 20 energetics of an animal (as calculated by the animal model 500 described below) to 
y population and community level phenomena (some of which may be included in the 
g calculation of other variables 800 step described below). For example, total activity time can 
be used to calculate certain life history variables like survivorship, growth and reproductive 
potential (because feeding frequency and handling time affect mass and energy intake by an 
25 animal, and thereby discretionary mass and energy), that are, in turn, used to calculate 
population dynamics. 

The temperature dependent activity/feeding frequency of an animal (i.e. the time 
during which an animal is active) is calculated as a linear function of the percentage of 
thermally available habitat (from the microclimate conditions output 400). A summation of a 
30 day's preferred activity times over a month and over the year then yields total annual activity 
time for the animal. 
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The total annual activity time of the animal is a key variable linking the individual 
energetics of an animal (as calculated by the animal model 500 described below) to 
population and community level phenomena (some of which may be included in the 
calculation of other variables 800 step described below). For example, total activity time can 
5 be used to calculate certain life history variables like survivorship, growth and reproductive 
potential (because feeding frequency and handling time affect mass and energy intake by an 
animal, and thereby discretionary mass and energy), that are, in turn, used to calculate 
population dynamics. 

The total annual activity time of the animal is calculated by summing the temperature- 
10 dependent activity/feeding frequency of the animal over all time intervals (over days, over 
months) totaling a year. The temperature-dependent activity /feeding frequency of the animal 
is determined by tracking the amount of time the animal is in its null condition as the animal 
model 500 is run - the null condition being initially defined by the user to be the normal 
activity /feeding time of the animal [e.g. whether the animal is active during the day (diurnal), 
15 night (nocturnal), dusk and dawn (crepuscular), etc.]. 

iv. Animal model conditions output 600 
The animal conditions output 600 may include many different types of calculations 
based on the user's preference. In general, the output 600 will include the values for several 
metabolic variables, the temperature-dependent activity/feeding frequency of the animal and 
20 the total annual activity time of the animal. 

The metabolic variables include the solutions generated using the animal model 500 
for several energy and mass metabolic variables (see FIG. 2). The heat balance solutions 
yield values for the metabolic heat generation needed to maintain the animal's core 
temperature (Qgen), the amount of heat entering the animal (Qin), the amount of heat leaving 
25 the animal (Qout). The mass balance equations yield values for the discretionary mass from 
food and water (Mdiscr) by calculating the mass of food and water absorbed by the animal 
(Mf,a, Mw,a), the mass of food and water ingested by the animal (Mf,i, Mw,i), the mass of 
food and water excreted and evaporated (Mf,d, Mw,u, Mw,evap), and the mass of food and 
water metabolized by the animal (Mf,met, Mw, met). Also included among the metabolic 
30 variables calculated is the lung air flow (including total air mass and C0 2 production). 

Optimal body size and other variables can be calculated for any given locality by 
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simply looping through the range of body sizes specified by the user and selecting the size 
that best meets user criteria, such as maximum discretionary energy (see System and Method 
of Implementation section below). 

d. Animal model conditions database 700 

5 The animal model conditions database 700 accumulates the animal conditions output 

600 over all time intervals, localities and species for which the models are run. The 
cumulative data are organized in the database 700 in such a way as to facilitate further 
calculations and/or display options. The location and time information continues to be 
associated with the animal conditions 600 as calculated by the animal model 500. This allows 

10 for the subsequent coupling of the animal conditions data set 600 with its associated temporal 
and GIS based information on climate, topography and vegetation. In this way, they can be 
used to generate spatially and temporally explicit interpretations of animal energetics and 
behavior, and displayed 900 or made available for further calculations 800 as described 
below. 

15 e. Calculation of other variables 800 

Further calculations may be performed using information in the animal model 
conditions database 700. When the animal model conditions data 700 are coupled with GIS 
based information on climate, topography and vegetation, they can be used to generate 
spatially explicit interpretations of animal energetics and behavior. 

20 For example optimal body size that maximizes discretionary mass or energy or other 

variables, optimal behaviors, intake of pollutants (e.g. pesticides), growth and reproduction 
potential as function of discretionary mass and energy survivorship/mortality as function of 
total annual activity time, and spatial variation in all variables (see for examples of these 
kinds of calculations, Adolph and Porter 1993, 1996, refs. above; Porter 2000, ref. above). 

25 The animal conditions data 700 may also be used to calculate community structure (if run 
over several species/body sizes). 

The analysis of spatial variations is also possible using the present invention. An 
example is given in the description of the system and method of implementation section 
below. 

30 f. Display of results 900 

The data in the animal model conditions database 700 and generated by any further 
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calculations 800 performed, may be organized and displayed in any number of ways. Any 
number of kinds of tables and graphs may be generated to display the data to best 
communicate the results of the model In particular, the data may be organized in tables 
particularly suited to use by spatial plotting software (such as ARCVIEW) to allow the data to 
5 be plotted in geographic space. 

g. Applications of the method. 

The method of the present invention can be applied to many and varied situations. 
Several examples follow. 

ii. Management and preservation of threatened or endangered species. 

10 One of the primary utilities of the present invention is to provide fully integrated 

microclimate and animal models which may be used in coordination to accurately predict the 
effects of climate and other environmental conditions on animals, that require relatively few, 
easily measured data to estimate food and water requirements, potential for activity time, 
growth, and reproduction for a wide variety of habitats. The fact that relatively few data are 

15 required is important to the present invention's application to animals in nature generally, 
because it eliminates the need for extensive sampling or smaller or otherwise difficult to 
sample populations. The present invention can generate the information necessary to identify 
potential reserves or transplantation sites, and to modify or manage existing habitats to 
optimize the growth, reproduction and survivorship of threatened or endangered species. 

20 Today, many species of animals are endangered. Often the cause of the endangerment 

is loss of their natural habitats to development or other land uses. Efforts to conserve the 
species require knowledge of the animal's adaptations to the environments in which they 
evolved, including their basic energetics. The question becomes one of how the animal will 
react both physiologically and behaviorally to the changes in their environments due to the 

25 changes in their natural habitats, displacement to new habitats, or other environmental 
influences. 

In this example, the integrated microclimate and animal models of the present 
invention were used to calculate the metabolic costs and food requirements of an endangered 
species of parrot, the Orange-bellied Parrot of Tasmania and Australia (Neophema 
30 chrysogaster, see, Porter et.al, 2000. Amer. Zoo., 40:597-630). With only approximately 200 
individuals still alive, reliable estimates of the metabolic cost, the necessary foraging time and 
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amount of food required to sustain the parrot in its seasonal habitats could lead to quantitative 
estimates of the amount of habitat needed to assure its survival. That information can then be 
used to guide conservation efforts to preserve or otherwise modify the environment of the 
parrot so that its discretionary energy and water, and activity time are optimized resulting in 
5 increased growth, reproduction and survivorship. 

Using the method of the present invention, parrot metabolism was calculated across a 
range of temperatures experienced by the bird. These calculations were then used to estimate 
the parrot's energy use requirements in its native habitat. The results for metabolic rate 
variation with changes in air temperature and with size of the bird (mass) are plotted in FIG. 

10 6. These results suggest that temperatures lower than 5 degrees C are stressful to hatchlings 
of lower mass. The method was also applied to calculate the effects of a range of foodstuffs 
ingested by the birds, on their growth and reproduction. These resulting data can be used to 
guide conservation efforts of the birds. 

ii. How vegetation alters the effect of climate on food requirements, 

1 5 geographic range, and predator-prey interactions. 

Vegetation can have an effect on the microclimates experienced by animals in addition 
to its role as a food source for some. Several of these effects are reported in Porter et al. 
(2001. Physiology on a landscape scale: plant-animal interactions. Unpublished manuscript). 
Below are summarized the reported effects of climate on food requirements and geographic 

20 range. Knowing these effects, managers of animals may modify environments accordingly to 
ensure the well-being of the animals. 

The present invention was employed to answer the question of how fires in 
Yellowstone National Park altered the nighttime thermal environment and energetic 
requirements of Elk overwintering in the park. Because the fires removed tree needles, the 

25 Elk were exposed to potentially cold stressful thermal conditions over the winter nights. The 
method of the present invention was employed to simulate the microclimates at the park in 
burned versus unburned locations, and the energetic costs and food requirements of similar 
sized elk compared. The results indicated that Elk in unburned locations required 20-25% 
less diurnal energy than did Elk in burned locations. The results further indicated that the 

30 differences in infrared radiation, not wind, resulting from the burn, was the major factor 
influencing the difference in energy requirements of the elk. Knowing these results, those 
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tending the Elk herds may be able to design feeding or other strategies to reduce the energetic 
stresses to the Elk, or otherwise modify their management practices. 

The way in which the thermal environment and vegetation (representing available 
resources) interact to determine the distribution of an animal, was investigated using the 
method of the present invention applied to the chuckwalla, Sauromalus obseseus, in the 
southwestern desert of the US and Mexico. The models of the present invention were used to 
calculate temperatures in the chuckwalla' s environments across the entire North American 
continent. Calculated patterns of degree days, coupled with GIS data on soil and vegetation 
types were compared to the known geographic distribution of chuckwallas. Plotting the 
results and comparing among the effects of temperature and vegetation/food, conclusions 
could be made on which factor(s) limit the animal's distribution. It was discovered that the 
northern limit of distributions may be constrained by egg incubation requirements and that the 
other boundaries of the distribution may be constrained by vegetation/food requirements. 

A plot of the degrees days available across the continent is shown in FIG. 7, with the 
chuckwalla' s range noted and the minimum degree days required for successful hatch of eggs 
specified (1705 degree days), by way of example. The plot represents the result of 
microclimate effects on soil temperature calculations for every degree latitude and longitude 
in North and Central America. Soil temperatures were calculated by the microclimate model 
of the present invention at the surface, and various depths below the surface. Only soil 
temperatures at a 15 cm depth were accumulated into a measure of "degree days." The 
calculations were for horizontal exposed soil with no shade on it. Successful chuckwalla egg 
incubation requires approximately 1700 degree days. Almost all of the current chuckwalla 
distribution lies in areas where the soil temperatures meet these minimum incubation 
requirements. If soil temperature alone is considered, the resulting plot suggests that the 
species should be able to live in Mexico and Central America. Why not? The other analyses 
suggest that vegetation type may be a critical variable limiting the chuckwallas distribution to 
the east, west and south. These results how the present invention may be employed to help 
manage the chuckwalla population in order to, for example, enable it to expand to the south, 
east or west, by modifying the vegetation/food available to it in those areas. 

iii. Other examples of applications of the present invention. 
Many other uses of the present invention are possible. Above are represented a few of 
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the uses relating to the conservation and management of animals. Those examples may be 
expanded to include applications to the management of biodiversity, zoo populations and the 
like. 

Other uses may include the analysis of disease vectors (e.g. mosquito vector) and how 

5 the effects of climate on their distributions may affect the spread of diseases. Still other uses 
include how pesticides may effect the energetics and survivorship of animals by quantifying 
the exposure of the animals; how the microclimate conditions experienced by animals stored 
in the cargo holds of airplanes are affected by the conditions they experience in order to 
increase survivorship by better controlling those environments (for example, carbon dioxide 

1 0 production by animals and oxygen availability as a function of the type of cargo hold, air flow 
and altitude); calculation of activity time, mass, and energy budgets for rare and endangered 
species to aid in their management (examples of species already investigated using the present 
invention include the giant panda of China, the Arabian Oryx on the Arabian Peninsula, the 
Yunnan snub nosed monkey on the China-Burma border, Puerto Rican parrots in Puerto Rico, 

1 5 marsupial gliders in Australia); assessment of the impact of manmade disturbances on the 
time, energy, and mass balances of rare and endangered species; forensic analyses to assess 
approximate time of killing of rare and endangered species in a natural habitat; environmental 
impacts of controlled burns in forests or national parks; assessment of climate change on 
distributions of animals in the past, present, or future; assessment of the impact of climate 

20 change on animal body size in the past present or the future; identification of regions of high 
natural selection pressures on animals in space and time; the design of animal shelters and 
animal husbandry practices in developing countries; and, dynamic changes in food web 
structure as influenced by climate, topography, and vegetation. This is a partial list of 
present and potential applications of the present invention. 

25 

DETAILED DESCRIPTION - METHOD & SYSTEM OF IMPLEMENTATION 

As described above, the method for calculating the effects of environmental conditions on 
animals of the present invention includes the provision of databases of environmental 
conditions 100 and animal characteristics data 190 for use in forming sets of environmental 
30 120 and animal 170 data to be input to the microclimate 300 and animal 500 models. These 
databases (100 and 190) may be very large (currently on the order of 9 terabytes/table for the 
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database MySQL). If used as input files instead of database tables, their size limits their 
accessibility by users who do not have access to larger mainframe computers. Such users are 
limited by the amount of memory and speed of their computers. The present invention using 
Practical Extraction and Report Language (PERL) and MySQL, addresses this problem by 
5 providing a method and system for its implementation on smaller computers, including 

personal computers. This allows many more users access to the methodology of the present 
invention and should enable its more widespread use among researchers, conservationists and 
other parties interested in understanding and managing animals. 

Referring to FIGS. 8 and 10, one embodiment of the system and method for 

10 implementing the present invention is summarized in a flow chart (FIG. 8) and basic system 
drawing (FIG. 9). FIG. 8 is presented to explain the basic steps in the method of 
implementation of the present invention that allow a user to access and process the large 
amounts of data necessary, using minimal RAM memory and disc space on a computer. 
Multiple modes of implementation of the present invention are, of course, possible. As 

15 depicted in FIGS. 8 and 10, the system is a client-server system, where the user's computer 
3000 is on the client side and the server 3010 may be located elsewhere on an Intranet or 
Internet. The same basic methodology could, however, be employed on a single computer 
with sufficient memory and disc space to store the databases and run the programs. 

The method comprises the generation of a graphical user interface (GUI) to the Web 

20 browser 3020 of the user's computer 3000 (1000 and 1010 in FIG. 8). The GUI can be 
customized to suit the particular application. For example, it can be customized for use by 
managers of endangered species to provide them choices and options that users in that field 
might typically want, or for use by environmental consultants. Likewise, an alternative GUI 
may be customized with features more appropriate to the researcher in a lab or field setting, or 

25 the like. 

In this embodiment, the GUI is provided in Hypertext Markup Language (HTML) for 
implementation with the client's Web browser 3020 and is generated by a PERL program 
running at the server. However other types of Common Gateway Interface (CGI) 
programming languages may be employed besides PERL, and the GUI provided in other 
30 languages besides HTML, The PERL driver, however, enables use of the databases rather 
than files for the inputs/ouputs, thereby minimizing the memory required and speeding the 
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processing. This, in turn, allows for the implementation on smaller computers. It also allows 
the use of massive sets of data, thereby enabling the spatial and high resolution temporal 
analyses possible using the present invention. 

The user then inputs a set of variables and/or conditions 150 to be used by the 
5 microclimate 300 and/or animal 500 models of the present invention 1020. A PERL program 
at the server (3030, representing all computer programs running at the server including the 
animal and microclimate model programs, output display programs, and the like) receives the 
user input 150 in HTML, or by direct manipulation of the PERL code 1030. The PERL 
program also reads one or more database tables contained in the environmental conditions 
10 database 100 and that supply input information about the environmental conditions (e.g., 
elevation, slope, aspect, vegetation type, and substrate type as a function of geographic 
coordinates) by locality 1040. After processing each locality input data set, the PERL 
program modifies the input data file to the microclimate model 1050. After each input file 
5 modification, the PERL program calls the microclimate program 300 to run 1060. Two 
yilS output files (called metout and soil) are created that form the above ground and below ground 
CO microclimate conditions 400 used as input for the animal model 500. As depicted in FIG. 8, 
Ifi the microclimate model 300 is programmed in FORTRAN, however, other programming 
u § languages may alternatively be employed. 

Q Depending upon how the user has specified animal data input (in the user-defined 

^ 20 variables and conditions 150), the PERL program will modify the animal input database 190 
^ to form the set of animal characteristics data 170 for input to the animal model program 500. 

The animal model program 500 is then run using the microclimate conditions 400 output from 
the microclimate model 300, and the animal characteristics data set 170, 1070. The animal 
conditions output 600 is stored in an animal conditions database 700, 1080, for later display 
25 900 and/or further calculations 800. 

If more time intervals remain, the animal model is run again 1070. At the end of all 
time interval runs, the program checks if the user has specified that other animal species are to 
be considered. If so, the program repeats from step 1040 until all species are considered. The 
program then likewise executes over all localities under consideration, by repeating step 1040 
30 on. 
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At the end of all runs, the data in the animal conditions database 700 may be accessed 
for further calculations 800 or for display 900. In FIG, 8, an example of displaying the 
results geographically is portrayed at 1090, in which the data are plot in GIS format using 
Arcview. These data are then provided to the GUI at the user/client 2000, and the plot, or 
5 other results, displayed. 

The user may specify that any number of calculations be made by the model at 1020. 
Depending on those specifications, the model(s) may repeat in several iterations, varying the 
data inputs each time. In this way a user may compare outputs for different body sizes or 
other animal variables to find, for example, an optimal body size or how variables vary 
10 spatially (described below), among many other possibilities. 

Optimal body size. Optimal body size and other variables can be calculated for any 
given locality by simply looping through the range of body sizes specified by the user and 
_ selecting the size that best meets user criteria, such as maximum discretionary energy. These 
sfj kinds of loops are in the driver code in PERL and are external to the FORTRAN executable 
1?= 1 5 elements called by the PERL code. 

Spatial analysis of animal energetics. The spatial analysis of animal energetics starts 
yQ by assembling data files that include digital elevation maps (DEMs). This process can often 
^ B be done by accessing files of quadrangle maps from various sources online such as various 
O United States government agencies. Alternatively, private corporations either have available 
Hj 20 or will create for a fee, digital elevation maps of anywhere on earth, thus enabling the input of 
^ geographic data for just about anywhere to the models of the present invention. Digital 
M elevation maps can be displayed in GIS programs, such as Arc View. Boundary maps or other 
geographic information can be superimposed on these images. ArcView can then translate 
this topographic image into an analysis of elevation. Slope and aspect (azimuth) of each 
25 "pixel" of information can be calculated as well These maps have a resolution of 30 m, so 
each "pixel" is 30 m on a side. These data can be exported into a text table. The collection of 
these data from separate tables into a single table that can be read by PERL (or other similar 
program) is easily automated. 

Each row of a table represents a unique location and the data associated with that 
30 location. For example, a row of data may contain latitude, longitude, elevation, slope, aspect, 
vegetation type, and soil type. Those data are then loaded into a database table. The PERL 
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program reads the input database table, one row at a time. Those data are used to modify the 
input data file for the microclimate model 300, basically modifying the geographic 
coordinates and other appropriate data in the input file to correspond to the particular location 
represented by the current row being read from the database table. Generally, a full year of 
5 microclimate simulations are run for each location and then the output is passed to the animal 
model 500. 

After all animal simulations are done for the same location, the PERL program reads 
another line from the database input table and repeats the process. After all calculations for a 
location are completed, the animal conditions output data 600 are written to the animal 
10 conditions database 700 output table. Its format is similar to the input data table, in that the 
first two entries in every row are the latitude and longitude followed by various animal output 
data for each month of the year and accumulated over the entire year. This output data table 
becomes the input table for Arc View or other GIS program that can plot discrete columns of 

C output as a function of latitude and longitude and display the information graphically. FIG. 9 

l7|15 is an example of such a plot showing the amount of discretionary energy (kJ/d) of a 60 kg 

Mountain Lion for the month of July, as a function of slope and aspect of the topography for a 

yO particular quadrangle of its range. Activity time, as well as other variables, may likewise be 

^ r " plotted spatially. 

fy 20 While the above description contains many specificities, these should not be construed 

!ff as limitations on the scope of the invention, but rather as an exemplification of the presently 
preferred embodiments thereof The above-described embodiments are set forth by way of 
example and are not for the purpose of limiting the present invention. It will be readily 
apparent to those skilled in the art that obvious modifications, derivations and variations can 
25 be made to the embodiments without departing from the scope of the invention. Accordingly, 
the scope of the invention should be determined not by the examples given, but by the 
appended claims and their legal equivalents. 
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